MD 4 Ca. SL AR sas. 


- 
ê 
: 
Ê 
- 
e 
= 


* 
4 Revista de Engenharia | 
ecnica Publicação da Associação dos Estudantes do Instituto Superior Técnico 


PROPRIEDADE, DA ASSOCIAÇÃO DOS ESTUDANTES DO INSTITUTO SUPERIOR TÉCNICO * DIRECTOR: M. DE ABREU FARO 
* COORDENADORES EDITORIAL; ISABEL ELIAS E HÉLIO SOARES * DIRECÇÃO E ADMINISTRAÇÃO: AEIST - Av. Rovisco Pais, 1000 
LISBOA - TELS.: 88 93 23/88 10 18 * PUBLICIDADE TÉCNICA - AEIST - Av. Rovisco Pais, 1000 LISBOA - TELS.: 88 10 18/89 75 15 
* COMPOSIÇÃO E IMPRESSÃO: CABOGRÁFICA, LDA., Rua das Biscoiteiras, Lote 18 (traseiras) - 2795 LINDA-A-VELHA * PUBLICAÇÃO 
TRIMESTRAL - 5000 EXEMPLARES * ASSINATURAS (4 números): PORTUGAL 500800 - ESTRANGEIRO 12 US$. 


técnica 3 « 88 NOVEMBRO 1988 


] 


tm 
EM 


EDITORIAL: 
ARTIGOS CONVIDADOS 


A ENGENHARIA DOS SISTEMAS DE DECISÃO 
L. Valadares Tavares 


ARTIGOS 


OS CRISTAIS LÍQUIDOS. SUAS PROPRIEDADES E APLICAÇÕES 
A. €. Ribeiro / J. L. Figueirinhas e M, H. Godinho 


A SÍNTESE DUM PROCESSO QUÍMICO INDUSTRIAL, UTILIZANDO A TECNOLOGIA DO PONTO 
DE ESTRANGULAMENTO 
Clemente Pedro Nunes e Bodo Linnhoft 


NOTAS CIENTIFICAS 


CONDUCTOR-BACKED DIELECTRIC SLAB: EFFECT OF UPPER-SHIELDING 
C. Paiva e A. M. Barboda 


POLYACETYLENE: A CONDUCTING POLYMER WITH ENORMOUS POTENTIAL APPLICATIONS 
M. A. G. Martins / A. P. B. Oliveira e C, A. €. Sequeira 


ISSN 0040-1714 


ERRATA DA REVISTA TÉCNICA Nº 2, PÁG. 59 e 61 


Expressão 
(2) 
(6) 


Lê-se 
— At2 ]ds=0 
À Bsd dl 


U! 


! 
Vo) 


Va 
Ki = -Ar/2 
Kgg =[ Jr? 
Ks3 = [-r By0--- 


r= m+(m2+ a2)1/2 


Deve lêr-se 
-— At]ds=0 


Elo. d B 


K1I1 =-A,º/2 
Ra DR 
Kog = [-+ BY9-— 


r= mt(m2+a2) 


EDITORIAL 


A saída a público deste número representa um indício claro de que a 
Técnica vai conseguir atingir a regularidade que se deseja e 
preconiza. 


De facto, ainda este ano sairá outro número e assim se conseguirá em 
1988 a publicação dos quatro números de 1988. 


Esperamos que em 1989 a Técnica ocorra com regularidade, trimestre 
a trimestre. 


E agradável registar que tanto bastou que a Técnica reaparecesse 
para que começasse, desde logo, a afluir mais original. 


Tenhamos pois esperança que nos serão asseguradas em 1989 as 
condições materiais de que actualmente dispomos. 


Tenhamos esperanças que se observará um aumento significativo de 
Notas Científicas e Artigos. 


lo ES i n dh SOCIEDADE DE REPRESENTAÇÕES, LDA. 


RUA DOS FANQUEIROS, 250, 2.º ESO. “| 875640/1 1100 LISBOA 


Não admira que já não se fabriquem carros como dantes. 
Fabricar uma caneta Já é bem difícil. 


Oh, que linhas elegantes do Hispano-Suiza de 1925! Oh, que linhas elegantes as da 


Parker Duofold de 1927! 


Há já muito tempo que o automóvel é inacessível, mas a caneta, regressa felizmente de 
novo. 

Desde há muito ansiávamos recriar este design prestigioso da Parker. 

O nosso próximo centenário forneceu-nos o desejado pretexto. 

Tal como os actuais “carros de sonho”,a Parker Duofold Centennial cinzela peças que 
são autênticas obras-primas. 

É por outro lado, orgulha-se de as produzir com tais cuidados que parecem fabricadas 
como nos velhos tempos. 

Em vez de moldar o corpo ea tampa, fizémo-los a partir de um bloco sólido, como então. 

Em vez de talhar o aparo em metal moderno, continuamos fiéis ao ouro. 

Em vez de fender o aparo com um qualquer dispositivo moderno, ainda o fazemos 
manualmente, utilizando uma lâmina mais fina que um cabelo. 

É um processo extenuante de fabricar uma caneta. Mas, tal como o Hispano-Suiza, o 


estilo e o prazer de usar são amplamente compensadores. 
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A ENGENHARIA DOS SISTEMAS DE 


DECISÃO 


L. VALADARES TAVARES 


IST. 1988 


RESUMO 


Apresentam-se neste artigo as principais linhas de 
desenvolvimento da Teoria da Decisão e suas 
contribuições para a Engenharia moderna. 
Discutem-se ainda alguns resultados inovadores 
obtidos recentemente no IST. 


ABSTRACT 


The purpose of this paper is presenting the major 
areas of development of the Theory of Decision 
Making and their contributions to the modern 
Engineering. Innovative results recently obtained at 
IST are also discussed herein. 


1 A ENGENHARIA ACTUAL: 
PARADIGMAS E EQUÍVOCOS. 


À semelhança do que acontece com a generalidade 
das noções importantes, o conceito de Engenharia 
não é fácilmente definivel. No que respeita ao seu 
objecto, talvez se possa dizer que as actividades de 
Engenharia visam controlar e aproveitar as 
condições e os recursos do meio ambiente, bem 
como realizar e gerir sistemas destinados à 
produção de bens ou serviços necessários às 
comunidades. Todavia, será provávelmente mais 
pelos paradigmas que propõe e adopta que a 
Engenharia se identifica e distingue doutras 
actividades. 


e nun... 


Original recebido para publicação em 27.04.88. 


A - Paradigma da racionalidade 

Em Engenharia, crê-se na obordagem pela Razão 
não se duvidando das vantagens que oferece em 
relação a outras abordagens, por exemplo, baseadas 
na afectividade. 


B - Paradigma da estabilidade 


Apesar de vivermos num mundo em permanente 
mudança e mutação, acredita-se na 
estacionaridade e, ou, estabilidade de muitos 
fenômenos o que permite identificar problemas a 
resolver e classificá-los segundo taxonomias diversas. 


C - Paradigma metodológico 


Julga-se prioritário estabelecer métodos para atacar 
as classes de problemas definidas desconfiando-se do 
livre arbítrio ou da simples intuição. Estes métodos 
assumem frequentemente a forma de sequência de 
estudos e regras ou até de manuais e regulamentos. 


D - Paradigma científico 


Crê-se nas potencialidades da Ciência para descrever 
o real em estudo, para compreender o seu 
comportamento e para melhorar a capacidade de 
intervenção. 


É - Paradigma da eficácia 


Apesar de não se ignorarem as múltiplas fontes de 
incerteza e aleatoriedade, acredita-se na capacidade 
de "engenhar" sistemas eficazes, isto é, sistemas que 
permitam atingir metas pré-definidas com índices de 
segurança ou fiabilidade muito elevados. 
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F - Paradigma da economia 


O desenvolvimento de qualquer actividade de 
Engenharia deve pautar-se por preocupações de 
economia no sentido de evitar qualquer dispêndio 
desnecessário de recursos. 


G - Paradigma tecnológico 


A realização de qualquer obra de Engenharia 
implica não só concebê-la e projectá-la mas também 
escolher e inventar materiais, dispositivos e 
processos que permitam concretizá-la da melhor 
forma. 


H - Paradigma comportamental 


A Engenharia afasta-se das Ciências Puras ou 
Aplicadas por não buscar propriamente teorias ou 
enunciados universais mas sim por pretender 
conhecer o comportamento da realidade em que vai 
intervir a fim de conhecer respostas em função de 
intervenções e cenários alternativos de aleatoriedade 
e incerteza. Esta preocupação "comportamentalista" 
reduz, evidentemente, a generalidade das suas leis 
acentuando a importância do caso especifico e 
atribui-lhe um pendor eminentemente pragmático. 


O desenvolvimento das actividades de Engenharia 
segundo estes paradigmas não tem deixado de criar 
equívocos diversos responsáveis frequentemente por 
atrasos e insuficiências no seu desenvolvimento. 


O primeiro e mais frequente diz respeito à sua 
relação com a Ciências que utiliza, em especial no 
que concerne à Matemática e à Física. É vulgar 
exagerar-se a sua identificação com ramos aplicados 
destas disciplinas esquecendo-se que os seus 
objectivos são sempre distintos pelo que a motivação, 
os estilos e a orientação não poderão ser 
semelhantes. 


Outro equívoco não menos habitual diz respeito à 
interacção com as Ciências Económicas que, 
tradicionalmente, se situam em áreas profissionais, € 
até escolas, distintas das da Engenharia. As 
preocupações de valoração e avaliação, de eficácia 
e eficiência já apontadas não são compativeis com 
atitudes “ separatistas " neste âmbito. 


Pelo contrário, a tendência actual tem sido a de 
procurar utilizar o contributo das diversas 
disciplinas pertencentes à área de Economia com 
intensidade não menor da relativa às outras Ciências 


Básicas. 


Por último, refira-se o equívoco da complexidade que 
tem vindo a favorecer abordagens mais complicadas 


em detrimento de outras mais simples só por aquelas 
aparentarem a utilização mais abundante de 
conhecimentos científicos. Na verdade, não é objecto 
próprio da Engenharia desenvolver teorias 
cientificas mas tão sómente aproveitar delas o que 
mais útil se pode extrair para melhorar as 
realizações da Engenharia. 


2. RACIONALIZAÇÃO DO PROCESSO 
DE DECISÃO 


2.1-0 MODELO DA TEORIA DOS SISTEMAS. 


Atendendo à já referida importância dos problemas 
de escolha ou decisão no contexto da Engenharia, não 
surpreende que se tenha vindo a estudar o modo de os 
formular, analisar é resolver evidenciando o muito 
que há de semelhante entre problemas de decisão 
vividos em especialidades distintas dos diversos 
ramos da Engenharia. Esta pesquisa tem-se baseado 
na Teoria dos Sistemas e na Investigação 
Operacional ( KAUFMANN, 1968) por oferecerem 
grau de generalidade apreciável e por permitirem 
representar muitas das facetas mais significativas 
(Ver, por exemplo, Bertalanffy, 1969). 


Os problemas de decisão a estudar prioritariamente 
apresentam fenomelogias muito variadas mas 
também algumas propriedades gerais a destacar: 


a) implicarem escolhas relativamente à 
afectação de recursos, humanos ou materiais; 


b) incluirem características e, ou, componentes 
não deterministicos cujo conhecimento é 
rodeado de aleatoriedade ou incerteza; 


c) apresentarem um carácter fortemente 
combinatório, ou seja, a configuração de cada 
decisão resultar da combinação de escolhas 
parcelares pelo que a extensão do conjunto de 
alternativas cresce rapidamente atingindo 
níveis muito superiores aos previstos pela 
intuição humana. 


O estudo destes problemas de decisão tem vindo a 
exigir um apreciável esforço de teorização, baseado, 
principalmente, na Teoria dos Sistemas a qual 
oferece a vantagem de permitir evidenciar o que há 
de essencial e geral na multiplicidade da 
diversidade,. Assim, importa recordar algumas 
noções assumidas como ponto de partida muito 
embora já resultem de operações de análise e de 
simplificação: 


Decisor (D) - Ente (singular ou colectivo) 
sujeito do processo decisório, ou seja, que goza 
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da propriedade de poder escolher, optar. Na 
verdade, este sujeito é frequentemente 
colectivo, designando-se então por actor todos 
aqueles que nele participam. 

Sistema (S) - Objecto da decisão, isto é. 
conjunto de realidades que serão afectadas 
pelas escolhas contidas na decisão o que 
pressupõe o estabelecimento duma fronteira 
g(S) que separe tais realidades das restantes. 


Alternativas (A) - Cada decisão possivel, isto 
e, Julgada merecedora de apreciação por ser 
realizável. 


“Acção (B) - Conjunto de actividades que 
concretizam a realização da escolha relativa a 
cada alternativa. 


Factores de incerteza (F) - Todas as 
solicitações que actuam sobre S e que não 
dependem das escolhas afectuadas por D. 


Resultado (R) - Conjunto dos efeitos que 
resultam, directa ou indirectamente, da acção 
inerente a cada alternativa. Designa-se por 
função de comportamento, (1), a relação 
existente entre R e (F;B). 


Informação do Decisor (1) - Todo o conjunto 
de conhecimentos e dados disponiveis pelo 
decisor nos momentos de decisão e que se 
supõem relevantes para as escolhas a efectuar. 


Modelo de decisão (M) - Representação 
simbólica de todas as operações relevantes 
para o processo de escolha em estudo. 


Estes conceitos permitem a imediata obtenção de 
certas noções complementares tais como a de espaço 
das alternativas (SA), espaços das acções (Sp) ou 
espaço dos resultados (0h) e bem assim a 
representação do processo decisório segundo o 
modelo mais simples da Teoria dos sistemas que se 
representa na Fig.l. 


2.2 O PROBLEMA DO GERAÇÃO DAS 
ALTERNATIVAS 


A definição de 24 não deve ser considerada como um 
dado a priori do processo decisório mas antes como o 
resultado dum processo construtivo que depende não 
só do espirito inventivo do decisor mas também da 
adopção de técnicas exploratórias. A pesquisa dos 
dominios de grandezas caracterizadoras das decisões 
e a combinação das opções distintas quanto às suas 
componentes têm sido aperfeiçoadas por técnicas 
diversas tais como a Análise Morfológica 
correntemente utilizada nas fases de concepção e 
inovação dos produtos de muitos ramos de 
Engenharia. 

Em geral estes processos e a própria realidade 
costumam originar espaços QA contínuos e com 
dominios muito vastos que não são fácilmente 
tratados racionalmente pelo que urge desenvolver 
técnicas de redução simplificadoras de OA . Esta 
redução traduz-se, em geral, pela discriminação de 
Qa substituindo-se por Q4' sendo Q4' um conjunto 
finito de alternativas representativas de QA. Esta 
representatividade pressupõe uma Taxonomia que 
sub-divida Qa.em classes, tanto quanto possivel 
representadas por um tipo único o que exigirá 
elevados índices de homogeneidade no interior de 
cada classe. Assim, Qa* será evidentemente, a 
colecção completa deste tipos. Observe-se que a noção 
de homogeneidade referida deve basear-se na função 
f já introduzida o que pressupõe conhecimento do 
decisor sobre S reforçando-se assim a perspectiva de 
que só em fases avançadas do processo decisório é que 
se vão obtendo definições aperfeiçoadas de DA ouD A. 


Em geral, mesmo utilizando as técnicas de redução 
referidas, é geralmente elevado o número de 
alternativas a estudar importando então 
representá-las de forma eficiente sendo de referir 
contributos disciplinares úteis tais como os 
produzidos pela Teoria de Grafos. 


O modelo apresentado permite evidenciar algumas 
dificuldades inerentes à análise dos conceitos 
apresentados: 


FIGURA 1 - Modelo do processo decisório 
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a) A distinção entre D e S não é fácil. 
Habitualmente desenvolvem-se relações de 
pertença ou rejeição com forte intensidade 
afectiva que não são compatíveis com esta 
distinção. Noutros casos, D é mesmo parte de S 
ou S parte de D. 


b) O conteúdo de D também não é simples. 
Para além de poder conter múltiplos actores, 
pode também ser considerado como uma 
estrutura de múltiplos decisores com 
alternativas e escolhas diversas. Entao, neste 
caso cada decisor deve incluir em S tudo o que 
respeita aos restantes decisores, 
designando-se estes problemas de jogo. 
Observe-se que não se estudarão estes 
problemas (RAIFFA, 1985) neste artigo por se 
crer que ultrapassa o seu âmbito. 


c) A distinção entre B e F também pode ser 
dificil de caracterizar. A conquista da 
capacidade de intervenção sobre S é, em geral, 
função também das escolhas de D. Outrossim, 
a capacidade de previsão por parte de De o 
próprio | reduzem a incerteza inerente a F, 


d) A noção de resultado implica distrinçar os 
efeitos da alternativa escolhida o que é, em 
geral, dificil, pois implica saber o que teria 
acontecido no caso contrário. À via resolutiva 
deste problema é a do estudo 
comportamental de S, ou seja, procurar 
conhecer a função R=T(F;B) referida, Op X Qp 
- Op. 


e) O fluxo de informação, |, resulta dos 
conhecimento de D, da redução de incerteza de 
F mas também do conhecimento de R que 
permite ir melhorando o conhecimento de f, 
Todavia, este melhoramento pressupõe uma 
noção de sequencialidade no processo de 
informação-decisão o que complexifica a 
formulação anterior. 


Todavia, os conceitos apresentados e o modelo 
representado são especialmente úteis por servirem 
de ponto de partida para a formulação das principais 
questões e linhas de orientação metodológica e que se 
enunciam nos sub-capitulos seguintes. 


2.3. O PROBLEMA DA TIPOLOGIA DAS 
DECISÕES 


Pretendendo-se racionalizar o processo de decisão a 
fim de apoiar o decisor, importa esclarecer que tipos 


8 


de escolhas existem já que têm consequências 
metodológicas distintas. Supondo que se determinou 
previamente QA', podem definir-se os seguintes tipos 
fundamentais de escolha: 

- Selecção dum sub-conjunto de Qa', Os: 
caracterizado pela sua preferência em relação ao 
seu complementar em Dar, Os. 


Esta noção de preferência aplicada aa! significa 
que qualquer elemento aeQs + é preferencial a 
qualquer outro b com be Qs-. 


As alternativas pertencentes a Qg: recebem 
designações diversas consoante o contexto Lais como 
"alternativas de mérito” ou “alternativas 
seleccionadas", enquanto que Qg é geralmente 
considerado como um conjunto de alternativas a 
rejeitar. 


- Escolha de uma e só uma alternativa 


Neste caso exige-se que Qs' tenha um só elemento o 
que, em geral, suscita mais dificuldades ao longo do 
processo decisório. Trata-se, pois, dum caso 
particular e mais restringido do problema anterior. 


- Distribuição dos elementos de Qa: por classes 
previamente definidas. 


Este tipo de problema corresponde à classificação de 
candidatos em classes de qualificação ou a problemas 
de afectação bem conhecidos onde cada classe 
corresponde a outras decisões bem definidas. 


Poder-se-à dizer que é uma generalização do primeiro 
problema que incluía apenas duas classes. 


- Ordenação dos elementos de Da". 


Este problema tem uma estrutura semelhante à do 
anterior exigindo-se que cada classe só contenha uma 
alternativa, ou, em casos excepcionais, permitindo-se 
que uma classe contenha mais do que uma 
alternativa desde que se considerem “empatadas”, 
isto é, equivalentes na apreciação que merecem pelo 
decisor. 


Como é óbvio, estes tipos de escolha podem estar 
associados a processos localizados num só instante do 
tempo ou, pelo contrário, a sequências decisórias em 
que os momentos de decisão são separados pela 
recepção de informação adicional, Podem também ser 
perspectivados para uma imagem (representação, 
conhecimento) única de S ou antes para conjuntos de 
imagens relativas e tudo aquilo que o decisor não 
controla, designando-se então a regra de escolha 
selecionada por D de estratégia. 
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2.4 O PROBLEMA DA CARACTERIZAÇÃO E 
QUALIFICAÇÃO DAS ALTERNATIVAS 


As alternativas devem ser caracterizadas por 
descrições que contenham, tanto quanto possível 
todas as perspectivas relevantes para a tomada de 
decisões. Assim, não interessará insistir em 
características não dependentes da escolha tais como 
aquelas comuns a todas as alternativas nem tão 
pouco acrescentar dados e informações resultantes 
apenas da preocupação de acumular pormenores 
sobre cada opção. Esta diferença de óptica afasta, 
frequentemente, o consultor do apoio à decisão em 
relação ao técnico especialista que se sente mais 
confortado com volumosas acumulações de descrições 
de especialidade. 


Entretanto, é geralmente ao longo do próprio 
processo decisório que se vão corrigindo e 
completando as caracterizações de acordo com as 
perspectivas de apreciação das alternativas que 
corresponde a eleger algumas características e, ou, a 
construir outras que se julguem importantes para a 
sua apreciação, as quais se designam de atributos 
ou de critérios ( no sentido mais geral ,ver ROY, 
1985). Estes critérios podem traduzir-se por 
grandezas definidas através de variáveis reais 
(investimento, caudal, distância, área, massa, 
centralidade, custo de manutenção, indice de 
poluição, frequência de produtos defeituosos, 
luminosidade, etc) ou por apreciações qualitativas 
dizendo respeito a perspectivas diversas tais como 
estética, tranquilidade, facilidade, simplicidade, 
harmonia, conforto, familiaridade, identidade etc.. 
Geralmente, é possível associar a cada uma destas 
perspectivas variáveis reais que contribuem 
quantitativamente para a sua definição já que se 
supõe que são causa da impressão produzida no 
decisor. Por exemplo, o conforto resultará também de 
certas variáveis climáticas, o belo de certos 
equilibrios formais, o simples de estruturas 
facilmente legiveis, a tranquilidade de ausência de 
perturbações, etc..Note-se, porêm, que este esforço de 
quantificação pode também ser nefasto quando 
conduzir a uma fragmentação exagerada de 
atributos, reduzindo o seu significado global. 


Muito embora sempre se tenha reconhecido a 
existência de critérios múltiplos na apreciação das 
decisões em Engenharia, a sua principal tradição 
consiste em representá-los por variáveis reais 
dimensionalmente homógeneas e assim compór uma 
função aditiva cujo máximo ou mínimo (sujeito a 
restrições) deve ser pesquisado, Esta abordagem tem 
potencialidades importantes (ver TAVARES, 1986, 
a) mas foi frequentemente adoptada de forma menos 
crítica ou esclarecida. Recentemente, 
desenvolveram-se metodologias particularmente 


úteis para o seu melhor aproveitamento (ver o 
sub-capitulo 3.1.) e outras que podem, 
alternativamente também ter grande utilidade no 
apoio aos processos decisórios (ver o sub-capitulo 3.2). 


25 O PROBLEMA DO ESTABELECIMENTO DE 
PREFERÊNCIAS 


A noção de escolha implica, óbviamente , a noção de 
preferência. Ora este conceito terá de ser aplicado 
segundo cada critério e terá de ser também 
estabelecido para a síntese dos diferentes critérios 
utilizados . Começando pelo caso mais simples, 
dever-se-ã permitir que o decisor possa comparar 
duas alternativas (a,b) segundo certo critério j e 
conclua, por exemplo, a P b, isto é, a é preferível a b, 
ou, pelo contrário, b P a. 


Isto significa que estabelecer relações de preferência 
é, fundamentalmente, introduzir o conceito de 
ORDEM em Q4'ºo qual se considera como uma noção 
primitiva (à semelhança das noções de elemento e de 
conjunto). 


Contudo, poderá também acontecer que o decisor 
considere estas duas alternativas como equivalentes, 
a É b, ou as considere indiferentes, a I b, muito 
embora não possa garantir a sua equivalência. Por 
último, poderá acontecer que o decisor não seja capaz 
de as comparar, a R b, ou que, pelo contrário, goste de 
graduar a intensidade de preferência, utilizando, pelo 
exemplo, dois ou mais graus de intensidade. Em 
resumo, compararmos realidades alternativas é uma 
atitude natural em todos nós mas não é de forma 
alguma evidente o modo de construir sistemas 
lógicos que traduzam tais comparações entre pares 
de alternativas. 


Estes sistemas dizem-se Sistemas Relacionais de 
Preferência (SRP) (ver ROY, 1985) e têm de tomar 
posição sobre diversas questões tais como o conjunto 
de relações binárias utilizáveis, as exigências de 
coerência, etc. 


Considerando que existe apenas P (e a sua inversa P, 
ou seja, se a P b então b P- a) tem-se o sistema 
chamado Ordem o qual é frequentemente utilizado 
em Psicologia (EDWARDS, 1954). Se se introduzir E 
diz-se Pré-Ordem ( VINCKE, 1980) e qualquer SRP 
permitindo a incomparabilidade diz-se parcial em 
oposição aos restantes que se podem qualificar de 
completos. A introdução de I origina novos SRP tais 
como o da Quase-Ordem ou da Pseudo- Ordem. 


Apresentam-se em 3.2 alguns dos seus conceitos mais 
importantes. 

À segunda questão consiste em saber como devem ser 
agregadas comparações binárias segundo critérios 
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múltiplos. É sem dúvida, uma questão mais dificil do 
que a primeira e que será abordada nos sub-capítulos 
seguintes. 


Todavia, importa começar por questionar o que 
poderá ser a racionalização desta etapa do processo 
decisório, julgando-se pouco controverso concluir que 
deverá assentar fundamentalmente na descrição do 
comportamento do decisor ( modelação de D). 
Todavia, se se reduzisse a esta preocupação, não se 
estaria a apolá-lo mas tão sómente a descrevê-lo, 
Ora, importa não esquecer que tais descrições ao 
serem conhecidas pelo decisor devem ajudá-lo a 
identificar-se, a esclarecer-se, a conhecer-se, 
estabelecendo-se para o efeito um processo 
interactivo Decisor & Modelo do Decisor. Na 
verdade, poderá definir-se um modelo interactivo 
como sendo todo aquele que vai alterando, 
interactivamente, o estado de conhecimento do 
decisor não apenas sobre S mas também sobre si 
próprio. 


Esta reflexão mostra, pois, que os instrumentos de 
apoio à decisão devem orientar-se não só para 0 
conhecimento de S mas também para o de D, o qual 
não deve ser encarado como um ente estático e 
pré-definido mas sim como um sujeito activo a viver 
um processo evolutivo. 


Esta nova orientação, distinta da tradicionalmente 
adoptada pela escola americana da Teoria da 
Utilidade ( FISHBURN, 1984), que pressupõe o 
decisor plenamente esclarecido sobre a sua estrutura 
de agregação de preferências não a alterando ao 
longo do processo decisório tem vindo a ser 
investigada principalmente no quadro das culturas 
europeias. 


Em resultado desta orientação, espera-se 
racionalizar não só o processo de escolha mas 
também o próprio decisor, considerando-se os 
instrumentos desenvolvidos como um apoio e não 
como um substituto do decisor. 


2.6 0 PROBLEMA DA VALORAÇÃO 


A sociedade moderna baseia-se intensivamente na 
representação de preferências por variáveis reais que 
se supõem quantificar a opinião do decisor. 
Considerando o caso unicritério, associa-se a cada 
par de alternativas (a,b), um real, v(a,b), positivo, 
nulo ou negativo, exprimindo a intensidade de 
preferência de a sobre b ( quanto mais elevado, tanto 
mais a é preferivel a b. Se fôr zero, a é equivalente a 
b e se for negativo, b é preferivel a a). Esta definição 
permite, obviamente, estabelecer relações de 
preferência ou equivalência entre acréscimos de 
preferências. ou seja, pode dizer, por exemplo que 


v(a,b)>v (c,d)>0 


o que significa que o melhoramento resultante de 
passar de b para a (b>a) é preferencial ao de (dc). 
Esta função v diz-se graduação. 


Observe-se que esta função é definida a menos de 
uma transformação linear e que, uma vez aceite a 
sua existência, é sempre possivel definir uma função 
real, valor de cada alternativa, v (a), também 
definida a menos de uma transformação linear e 
sendo 


v(a,b)=V(a)-V(b). 


Esta quantificação corresponde a traduzir a 
intensidade de preferência de a em relação a b 
através da distância euclidiana definida em Ri entre 
V(a) e V(b). Tem-se pois, a habitual condição de 
cadeia: 


d(a,b) = d(a,c) + d(cicg) +... + d(cm,b) 


representando por d(x,y) a referida distância entre x 
eyeporta, cj cz, ..Cm,b) um conjunto de alternativas 
em estudo. 


Observe-se que a escolha da métrica euclidiana é 
plenamente arbitrária talvez resultante da sua 
utilidade no estudo do mundo real onde vivemos. E, 
todavia, mesmo adoptando esta métrica, convirá não 
esquecer que a condição de cadeia apresentada já não 
é válida em RN associando as N dimensões deste 
espaço aos N critérios de apreciação pois então, em 
RN,d(a,b)sd(a,c;))+ d((cy co) + d(cycg) +... +d(cn,b) 


É muito duvidosa a justificação desta inequação se 
aplicada aos problemas de preferência pelo que os 
defensores das métricas euclidianas aplicadas a estes 
problemas preferem ficar em Rli(ver, por exemplo, 
HICKS et al, 1943). 


A fundamentação associada à existência desta 
função valor para a decisor individual tem vindo a 
ser estudada especialmente na Teoria Economica 
desde o século XVIII (BENTHAM, 1789). Na 
verdade, numerosos autores têm considerado que a 
escolha individual está associada à maximização 
duma função definidora da sua satisfação em 
resultado da alternativa escolhida supondo-se que o 
“homem-económico” possui informação completa, é 
infinitamente sensivel e adopta uma abordagem 
racional (Neste contexto, V (.) designa-se por função 
utilidade ). 
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Infelizmente, a aplicação deste conceito à 
combinação de diferentes bens levanta dificuldades 
sérias, surgindo então a hipótese da aditividade 
linear para facilitar a sua resolução ( EDWARDS, 
1954). Todavia, na generalidade dos casos, esta 
hipótese é rejeitada por fortes evidências 


experimentais o que justificou o conceito de curva da 


indiferença (EDGEWORTH, 1881) que não é mais 
do que a partição de Qa' em 9,,....Qm tendo-se 
Q,N9;=0 com izj (i,j=1,..., M), U;Q;= Q4:, supondo 
que se pode sempre estabelecer uma relação de 
ordem em Q,,.....Qm o que significa que Q,é 
preferivel a ou é preferido por O; com izj,i,j=1,...., 
M. (Note-se que não é permitida a equivalência entre 
Q, e 92;). Esta evolução corresponde a deixar de ter 
uma definição cardinal da função utilidade passando 
a adoptar-se uma definição baseada na relação de 
Ordem já referida. Observe-se também que a 
generalização de V (definido em R1) para o caso do 
decisor incluir grupos sociais tem sido estudada na 
Economia do Bem Estar (GRAFF, 1957) tendo-se 
inicialmente procurado estabelecer uma função 
utilidade aditiva (KALDOR, 1939) que se veio 
mostrar fonte de equívocos e incoerências. A 
moderna Teoria do Bem Estar (MISHAM, 1981) tem 
vindo a explorar as potencialidade de SRPs sendo 
importante abservar que o Teorema de 
Impossibilidade de ARROW (ARROW, 1951) vem 
sendo, recentemente, resolvido através da 
flexibilização das condições de coerência do SRP 
adoptado. 


Perguntar-se-à então o que se pode concluir quanto à 
existência eventual duma função valor se se adoptar 
uma abordagem baseada em SRPs. A resposta a esta 
questão tem sido obtida durante os últimos 20 ou 30 
anos, podendo dizer-se que em certos casos (por 
exemplo, a Pseudo-Ordem, a Quase Ordem e a 
Hiper-Ordem) existem funções reais, g(a), 
interpretáveis em termos de intervalos, o que 
significa, por exemplo, (g(a)-g(b))>p equivale a um 
determinado juizo preferencial (VINCKE, 1980). 


2.7 O PROBLEMA DA AGREGAÇÃO 
MULTI-ATRIBUTO 


Como já se referiu, a generalidade dos problemas de 
decisão implica a consideração de múltiplas 
perspectivas para qualificar as alternativas 
constantes de 924". Estas perspectivas não deverão ser 
consideradas como "primitivas" na racionalização do 
processo decisório mas sim como o resultado de 
operações prévias de abstracção e análise pois cada 
atributo é o resultado de definições mais ou menos 
elaboradas (fiabilidade, custo, segurança, conforto, 
dedicação, eficiência, economia, pontualidade, etc.). 


O número de atributos depende muito do problema 
em estudo mas julga-se que muitos dos problemas de 
Egenharia podem ser estruturados com base em três 
atributos, (TAVARES, 1984) cada um dos quais pode, 
por sua vez, resultar de outros sub-atributos. 
Citem-se, a título exemplificativo, algumas 
situações-tipo: 


a) avaliação e selecção de empreendimentos: 


. Custos de investimento 

. Custos de funcionamento e manutencão 

. Serviços ou outros benefícios gerados 
pelo empreendimento 


b) escolhas tecnológicas: 


. Custo-eficácia (viabilidade de curto prazo) 
. Impactos ambientais 
. Efeitos sociais e politicos 


c) Gestão de projectos 


. Custo-beneficio actualizado 

. Prazo 

Risco de não cumprimento das metas de 
calendário 


Para certos atributos, não é dificil associar-lhes uma 
função valor, V;(.), enquanto que para outros apenas 
se podem estabelecer comparações entre pares 
alternativos. 


A forma de agregar as comparações e, ou, valorações 
segundo atributos múltiplos tem sido estudada por 
numerosos autores propondo modelos, métodos e 
técnicas diversas, (ver por exemplo, ROY 1973 e 
SAATY, 1977). Em termos gerais, podem 
classificar-se estes processos em três grandes grupos: 
os compensatórios, os não compensatórios e os mistos. 
No primeiro grupo, admite-se sempre a possibilidade 
de estabelecer compensações entre os diversos 
atributos sendo nomeadamente a média ponderada 
das funções de valor para os M atributos o exemplo 
mais conhecido e adoptado. 


Pelo contrário, a abordagem não compensatória exige 
relações ou qualificações especificas para cada 
atributo ( e independentemente do que se verifica nos 
restantes atributos) a fim de poder concluir sobre o 
valor ou a preferência entre alternativas. 


E., por exemplo, o caso de num concurso com N 
candidatos, e M atributos, atribuir a cada candidato 


uma classificação, V (i), dada por 


Vh)= min ; V; (1) ou V()= max ; V; (1) 
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Semelhantemente, poderão estabelecerem-se 
comparações binárias entre pares de candidatos (por 
exemplo, a sobre b) atribuindo a comparação menos 
(ou mais) vantajosa para a segundo os M atributos 
considerados. 


Os métodos mistos incluem análises dos dois tipos 
anteriores e não se duvida ser o modelo Electre | (ver 
ROY, 1973) aquele que mais se tem divulgado. Neste 
modelo, diz-se que a prevalece ' sobre b (a não é 
preferivel em relação a b) se e só se o indice de 
concordância de a sobre b (C,p) fôr superior ou igual 
a certo limiar (1, limiar de concordância) e ainda se o 
indice de discordância (D,p) fôr inferior a certo 
limiar (1 limiar de discordância). 


O índice de concordância é de 
compensatória já que é dado por 


natureza 


sendo J o sub-conjunto de atributos para os quais b 
não é preferivel a a e sendo W; o coeficiente de 
importância relativa do atributo j. 


Pelo contrário, o indice de discordância é não 
compensatório, assumindo o carácter de veto: 

Da = max | Vj(a)- Vj(b) | 
je dJ' 


sendo J'o sub-conjunto de atributos para os quais 


Vila) < V;(b). 


2.8 O PROBLEMA DA INCERTEZA E DO RISCO 


Um dos desafios mais polémicos e sempre lançados à 
racionalização do processo decisório consiste no 
tratamento da incerteza e do risco. Na verdade, a 
existência do vector F e a dificuldade de predizer o 
comportamento de S originam dúvidas sobre R pelo 
que cada elemento de Qa* passa a ser caracterizado 
por mais do que uma configuração alternativa de 
atributos, agora designados de predicados atendendo 
à sua variabilidade. Em certas situações, ditas de 
Risco , é possivel probabilizar as realizações 
alternativas das configurações referidas podendo 


(1) O termo prevalência é proposto pelo autor para 
traduzir o termo inglês outranking ou o equivalente 
francês surclassement. 


então tirar-se partido da abordagem estatistica 
(Teoria Estatistica da Decisão, Raiffa, 1970). A 
introdução desta hipótese na valoração das 
alternativas foi especialmente estudada no quadro 
da Investigação Operacional devendo referir-se a 
conhecida Teoria da Utilidade de Von Newmann e 
Morgenstern baseada nos axiomas que se 
apresentam seguidamente (Von Newmann et al, 
1954). 


Designando-se por L a participação do dominio de V, 
2, num conjunto L=1,,...... lesendo |;Nl;=0 com iz) 
e U;|;=9,, e designando-se por pa(i) a probabilidade 
de ocorrer |; tomando-se a decisão a, admite-se que: 


Axioma 1: Para qualquer par de alternativas 
(a,b) caracterizados por pali) e pp(i), tem-se 
sempre um dos três casos. 


aEb(ebEa)ouaPbou bPa. 


Axioma 2: E e P são transitivas, ou seja, se 
abbebEc(ouaPbebPc)entãoaE c(ou 
aPc). 


Axioma 3 :Para qualquer trio de alternativas 
(a,b,c)comecPbeb Pa; existirá um valor 
único de IL ( O=Il<1) para a qual a lotaria 
((a, e); II) é equivalente a b. (Designa-se por 
lotaria " (a,c); 1 " uma nova alternativa tal 
que, se fôr escolhida, corresponda à realização 
dum ensaio aleatório do qual resulta a 
ocorrência do resultado a com probabilidade Il 
e do de e com probabilidade (1 - 11 )). 


A inclusão nesta abordagem da hipótese 
multi-critério tem sido também objecto de trabalhos 
recentes ( ver FISHBURN, 1977). 

Outra linha importante também desenvolvida para 
tratar os problemas de risco diz respeito à Teoria dos 
Extremos que concentra a sua atenção no estudo dos 
eventos raros (baixas probabilidades) mas altamente 
significativos pelas suas consequências ( catástrofes, 
grandes sucessos, etc.). 


Estas duas perspectivas têm vindo a tirar partido da 
noção de esperança matemática cujas 
potencialidades e limitações foram inicialmente 
discutida por Bernoulli. 


3. DESENVOLVIMENTOS RECENTES 


A fim de ilustrar as potencialidades das metodologias 
de apoio à decisão usando as linhas de racionalização 
apontadas anteriormente, é interessante apresentar 
dois resultados importantes recentemente obtidos no 
Instituto Superior Técnico. 
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3.10 MODELO TRIDENT (in TAVARES, 1984) 


Conforme se referiu anteriormente, a adopção da 
utilidade compensatória na agregação, para cada 
alternativa i, das suas apreciações segundo os 
diversos critérios ( j=1, ....,M), tem-se traduzido 
geralmente pela definição duma função utilidade 
(dita " Função Utilidade Multi-Atributo", in Roy 
1985) que agregue as funções de valor relativas aos 
diferentes atributos ( V; (1)) : 


V dy) — U(V | + à Pa V; e V M (1),). 


Ora, o caso mais frequentemente adoptado consiste 
em admitir para U uma estrutura linear e aditiva, ou 
seja: 


V (1) = 2W, uu, (1) 
Sendo Wjos coeficientes de ponderação adoptados 
(com L; uj=1 ) e representado por uj(i) uma 


transformada de V; (i). Esta transformada pode 
assumir formas diversas tais como : 


a) uj(i) = Vj(); 


VD) — min, (V (0) 


6) us) =: ———eeeeeeeeee 
É max V(D-min. Vi) 
e LJ 


se se pretender padronizar a contribuição de cada 
função valor em termos relativos da variação 
conhecida em Dar; 


ce) u; (1) = f;(V;6D)-min; V; ()) 
ou 
c)u;(i) =-f; (max; V;6)-V;G)) 


se se pretenderem também arbitrar " razões de troca” 
("trade-offs") para inicio de análise: 

Em 

m m 

feW, 
A escolha da tranformação mais aconselhável 
depende, como é evidente, da natureza e das 
características de problema em estudo. Todavia, 
surge sempre como questão central a escolha da 
ponderação pelo que importa desenvolver modelos 
que permitam aprovar a escolha dos coeficientes W. 
Assim sendo, e atendendo à importância especial de 
que se reveste o caso de M ser igual a 3, 
desenvolveu-se um modelo original (“Triangular 
Decision Technique" - TRIDENT) aproveitando 
também conceitos da abordagem baseada em 
comparações binárias, que, e através dum processo 
de interacção com o decisor, o ajuda a tomar posição 
quanto à ponderação. 


ôu (DJ yu =" 


Este modelo, após algumas análises preliminares, 
estuda o problema da ordenação das N alternativas 
(i=1,...,N) determinando no espaço dos pesos 
(definivel em R2 já que W/+ Wo + W3=1)os lugares 
geométricos da equivalência e da preferência para 
cada par (i, ii). 


W2 


0 | W 
FIGURA 2 - Espaço de ponderação 


] 


O espaço das ponderações possiveis (1w) é limitado 
pelo triângulo apresentado na Fig 2ea 
equivalência entre duas alternativas (1, 1') é 
representável por uma recta ( recta de equivalência) 
atendendo à expressão analitica de V (1), obtendo-se 
para os N ( N-1)/2 pares possiveis "mosaicos" do tipo 
do incluido na Fig. 3. Observe-se que cada recta de 
equivalência, se interceptar Q w, particiona-o em 
dois sub-espaços (Q'w, 2'w) sendo uma das 
alternativas preferível à segunda num deles e 
preferida pela segunda no outro sub-espaço. Assim se 


0 I W, 


FIGURA 3 - Mosaico das rectas de indiferença 
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conclui fácilmente que cada sub-espaço não 
interceptado por rectas de equivalência ((9k)) é o 
dominio de validade de cada ordenação possivel para 
Qapelo que a sua localização e a sua extensão dão ao 
decisor informações que poderão ser muito úteis para 
concluir sobre a adequação e a estabilidade de cada 
ordenação. Estes sub-espaços Qp, são 
necessáramente convexos e a sua caracterização e 
determinação podem ser fácilmente implementados 
computacionalmente 


QUADRO 


Dessto | [ic [Ro 
(=p 


CARATERIZAÇÃO DAS 
ALTERNATIVAS POSSIVEIS 


À titulo exemplificativo cite-se o problema de 
planeamento de certo empreendimento (TAVARES, 
1986) para o qual se constroem quatro decisões 
alternativas caracterizadas pela duração total sem 
ocorrência de atrazos, T (em semanas), pelo valor 
actualizado do custo total ( em 105 dolares dos EUA), 
e pela probabilidade de não conseguir completar o 
empreendimento de T devido a atrazos, risco R (em 
percentagem). 

À presentam-se no Quadro 1 as características destas 
alternativas. Adoptando-se a transformação c) com 
fi=f,=f9=1, obteve-se o Quadro 2 e o mosaico 
apresentado na Fig.4,pormenorizando-se no Quadro 3 
as razões de troca para os vertices dos poliedros Ox. 


QUADRO 2 


CORRE 


MARIZ u j(i) 


Em função das razões de troca referidas. o decisor 
considera apenas aceitável o sub-espaço (a,b,c,m) a 
fim de evitar qualquer ponto com 


du (0) du (1) 


vais OU Vi) 


Bu li) Bu li) 


igual ou superior a 2. Ora, com esta restrição, 
verifica-se que as alternativas vencedoras são 1 ou 2, 
preferindo-se 2 no caso de W3 ser menor do que 71% 
do peso atribuido ao segundo critério. Esta 
consideração colhe a adesão do decisor que acaba por 
seleccionar 2 deixando 1 em segundo lugar. 


Aplicações deste modelo a outros dominios 
(designadamente, da Engenharia Municipal). 


QUADRO 3 
ANÁLISE DE SENSIBILIDADE 
Ponto (óu,(4)/%u 


I - 
2 0 o 
3 0 - 
“ - 0.67 -— 
5 - 1.05 — 
6 - 2.02 - 
7 - - 
- - 13,84 
9 - - 1.31 
IO - - 0.71 
1 - - 0,48 
12 - 1.05 - O.71 
13 - 1.05 - 1.38 
14 - 2.02 - 1.31 
15 - 2.02 - 0,73 
16 - 1,38 - 0.93 


(Suy(i)/%u,(1)) em sem./10º u.s.dol. 


(duy(1)/6u, (1) em x p10º u.s.dol. 


encontram-se já publicados em revista da 
especialidade ( ver, por exemplo, BANA e COSTA, 
1988). 


Observe-se, portanto, que este modelo é , no fundo, 
uma análise de sensibilidade desenvolvida no espaço 
dos pesos e será concretizada computacionalmente 
oferecendo diversas iteracções gráficas e numéricas 
como o decisor a fim de o ir ajudando a conhecer 
melhor as alternativas a as suas próprias 
preferências. 
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ORDENAÇÃO OPTIMA 


D.056 
0.016 
0.041 
0.06 

0040 
0027 


0008 
0028 
0023 


E) 
ao Split) Vul 


0003 
0004 
0004 
0.003 


FIGURA 4 - Mosaico obtido pelo modelo TRIDENT 
para o problema do quadro 2 


3.1 A HIPER - ORDEM (IN TAVARES, 1988). 


O SRP tradicionalmente mais utilizado é a 
Pré-Ordem que inclui, como já foi referido, as 
relações de equivalência e de preferência. Todavia, 
este sistema apresenta importantes limitações, 
especialmente resultantes da transitividade de E tal 
como foi claramente explicitado pelo famoso 
paradoxo, de Luce ( ver Luce, 1956): Em resumo, 
tendo-se a E b;, by E by,....b, E c, é fácil encontrar 
situações em que a E c não é verdadeiro. 


Por estes motivos, surge a relação de indiferença. 1, 
que já não é transitiva, podendo dizer-se que o mais 
conhecido SRP desenvolvido e incluindo | é a 
Pseudo-Ordem (Roy, 1973) que inclui também as 
relações de preferência fraca (Q) e forte(P) 


Supondo que se reduzem todos os conjuntos de 
alternativas equivalentes entre si a uma única 
decisão pode-se simplificar a apresentação deste 
sistema omitindo E. Tal como é habitual, 
considera-se aqui que | é simétrica e reflexiva, Qe P 
são assimétricas, e, portanto não reflexivas. 


Neste sistema, as condições de cadeia exigem que: 
SIS. Cs 

S2nI2=0 

PRPECP 

P2n P2=0 
PILIQCP,QIPPCPQICP,LQ.P.CP 


utilizando-se as convenções seguintes ( sendo R, 
Ry,R9,R3,R4 elementos ou sub-conjuntos do conjunto 


(P-, Q-. 1. Q. P)): 


a) R;.Ro.(ou Ry,R9,R9) CR4 
significa que a Rb com RCRg se (aR,,b) N (bR9c) (ou 
(aR; b)N (bRoc) N(cRgd )); 


b) aR;2b significa (aR;c N cR;b); 
c)JS=QUP 


dja R b significa que não é verdade, quer a R 
b, quera Rb, quer bR a. 


Este sistema apresenta potencialidades 
interessantes, especialmente por se poder 
demonstrar que existe sempre, para dado conjunto 
finito de alternativas,  Qa: uma função 
Pseudo-Critério, g(a) sendo a o elemento genérico de 
Qa'e tal que (VINCKE,1980): 


alb « -q <g(a) <q 
aQbe q <gta) - g(b) <p(g(b)) 


aPb <p (g(b)) Zgt(a) - g(b) 


em que q (não negativo) é o limiar da indiferença e 
p(g(.)) é uma função real definida sobre Qa: 
satisfazendo: 


p(gt)Jizq 
p(g(a))-plg(b)) 
gta) - g(b) 


Esta função é uma generalização da noção de função 
de valor já apresentada sendo, porêm, o significado da 
diferença | g(a) - g(b)] dado em função dos escalões q e 
p, ao contrário do que acontece com a função valor 
ou utilidade. Estes escalões podem ser considerados 
como resultantes de limiares de sensibilidade do 
decisor, abandonando-se, pois, a hipótese da 
sensibilidade infinita. 


Infelizmente, as condições de cadeia apresentadas 
são dificilmente justificáveis e alguns dos resultados 
que permitem não serão fácilmente aceites pela 
generalidade dos decisores tais como (alclb)jN(a 
P b ) mas, em contrapartida, excluem outros 
possivelmente compativeis coom a lógica de muitos 
decisores: 
(aPbPejN(aQbQe) 


(aPbQeQd)Nlald) 
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Em conclusão, julgou-se importante construir um 
outro SRP também incluindo as relações (P- QI, Q, 
P ) mas mais flexivel do que a Pseudo - Ordem por 
possuir condições de cadeia mais simples e mais 
gerais. Com este fim, generalizou-se a noção de 
transitividade, introduzindo-se o conceito de 
Transitividade generalizada (TAVARES, 1987) que 
consiste em poder exprimir todas as condições de 
cadeia por cadeias com um número mínimo de 
comparações, ou seja, à R; b, bR9c. como é evidente, 
a Quase - Ordem não possui esta propriedade. 


O sistema que se propõe ( Hiper - Ordem ) possui 
esta propriedade e as suas condições de cadeia 
resultam dum confronto de condições e principios de 
racionalidade considerados "primitivos". As 
condições lógicas são as seguintes : 


| - Comutatividade 


Se Ri.R9=Rg e 2 ( Rj, R9) (representando-se 
por (R4, R9) o sub-conjunto de P. Q,I,Q, Pao 
qual Rg deve pertencer), e se R$ R;=Rg€e Q' 
(Ro, Rj) então deverá ter-se = 9" 


2 - Inversão 


Se R,.R5=R5 e Q ( Rj, R9) ter-se-à 
Ri=R9=R4= e OQ (Ry, R9) em que Q é o 
conjunto das relações inversas contidas em 2. 


3 - Coerência 


Deve-se garantir a existência de Ry «e O (Ri, 
Ro) para a qual Ryce O (R5, Ro)ede Ryce O(Ry, 
Ro) para a qual R$€e O(Ry,R9) 


Os princípios estabelecidos para obter as condições 
de cadeia são as seguintes : 


|- Princípio do Reforço 


Se Rj e R9 representarem preferências do 
mesmo sentido, Ry=R;,.R9 também o deverá 
ter com intensicade não inferior à da mais 
intensa do par (R,, R9): 


P.P -»P 
P.Q>P 
Q.Q =>Q:P 


2 - Princípio da Indiferença Duvidosa 


Se uma das relções (Ry, Rg) fôr I (por exemplo, 
R,), então Rg=R,.Rg poderá ser Rg ou uma das 


relações que lhe fôr adjacente no conjunto 
ordenado P:, Q,1,Q,P: 


PI>Q: P 
QI>P:QII 
L1->Q;1;Q 


3.Principio da Oposição 


Se Rj e R5 tiverem sentidos opostos e igual 
intensidade, R3 poderá ser R;,, Ro ou 
umarelação compreendida entre Rj e R$ no 
conjunto P,Q,1,Q,P: 

PP >P-: Q:.1.Q:P 

Q.Q ->Q:1,Q 


4 - Princípio da Oposição Atenuada 


Se R; e Ry forem preferenciais de sentido 
oposto e intensidade diferente, Ryg=R,.Rs5 
estará compreendida entre R, e R5 no conjunto 
Q ou será a mais forte do par (Ry, R9): 


PQSLQP 


Em consequência destes principios, obtêm-se as 
condições de cadeia seguintes : 


laPbPe>aPe 

2)JaPbQcetlou aQbPc)-»aPc 

3)aPbl c(ou albPc)->»aQc; aPc 
4jaPbQ- clou aQ-bPc)->alc; aQe; aPc 
5)aPbP-c(ou aP-bPc)->aP-c; aQ-c; alc; aQc; aPc 
6)aQbQe—aPc; aQe 

TaQble (ou albQc)->alc; aQe; aPc 
8)JaQbQ-clou aQ bQc)->alc; aQ-c; aPc 
9)aQbP-clou aP-bQc)-+aP-c; aQc; alc 
I0)alble->aQ-e; ale; aQe 
IDalbQ-e(ou aQ-ble)->aP-c; aQ-c; ale 
12)albP-elou aP-ble)->aPc; aQ-ec 
13)aQbQc>aP-c; aQe 
14)aQ-bP-clou aP bQ-c)>aP-c 
15)aP-bP-c>aP-c 


É agora necessário verificar se as condições lógicas 
são satisfeitas por estas condições de cadeia. Esta 
comparação é imediata para 1 e 2 podendo 
verificar-se que também a condição de coerência é 
sempre satisfeita através do Quadro 4. 


Este sistema apresenta potencialidades valiosas, 
designadamente as seguintes: 


a) utilizando alguns teoremas da Teoria dos 
Grafos (ROY, 1970) pode demonstrar-se a 
possibilidade de representar numéricamente 
qualquer comparação, utilizando agora uma 
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QUADRO 4 


VERIFICAÇÃO DE COERÊNCIA 


Caso R | k> Rs = (RR) MR con k ha R, E (8) com R oo dh, Verificação 
segundo as coníiiçoes ne Op) e O(R') 
de cadeia 
| P p P.Q,I,Q.r P",Q.1,Q,P Sin Sia 
2 Q p I,Q,P PO .1,0,P Sim Sim 
1 | Q ou ! 1,0.P, ou O,P PQ .I ou PL.Q.LO! Sim Sim 
“ P Q. I ou QouP [,0.P ou Q,P vu P Q, P” ou PQ LI, ou 
FP O LAP 5 im Sim 
5 P P P'ouQ oulouQouP P",Q,L,Q,P ou P” ou P ou Q,M” ou Sim Sim 
1,Q,P, ou Q,P ou P',Q, ou 
P ou P PQ LAP 
6 q Q Q ou P Q",.1,Q ou 1,Q,P Q",1,0 qu 1,4,P Si Sim 
DS O SG = Soco sons om neea pie io op senao ia ooo ra eo dias fecano se 
7 Q | ou Qou P Q.1,0 ou T.O,P' PT ou Q,I,0 Sim 5 im 
ou Q,P ou 1,0, 
H Q o Q cuTouQ Q,L,Q ou LOP PQ ou PQ, “Sim Sim 
ou Q,P ou 0 ,1,0 
LO I | Gp uu Tou Q PATA ou VAO PO A ou PIA Sim Sir 
ou [,Q,P ou 1,0, 


CSS 


nova função, dita função de comparação f(a,b), 
associada a cada par de alternativas, Esta 
função é também interpretável em função de 
limiares de indiferença ou preferência (q ou p) 
as quais agora podem ter de variar de par para 
par de alternativas. Na verdade, esta função é 
uma generalização da função Pseudo - 
Critério, distinguindo-se desta pela sua 
definição sobre (24'XQ4') em lugar de o ser 
sobre Qao que é claramente mais coerente 
com a orientação adoptada de basear a analise 
de Qa' nas comparações binárias dos seus 
elementos. 


bja variabilidade de (q,p) permite introduzir 
heterogeneidades no grau de sensibilidade do 
decisor passando agora a ser possivel ter ( 
boQbr; byQbo ,...., baiQb )NbolcNIb,) tal 
como se representa na Fig. 5. Esta situação 
pode ser exemplificada por casos que 
generalizam o paradoxo de Luce ( ver LUCE, 
1956 e o paradoxo de TAVARES, 1987) e que 
assumiam a forma de incoerências no quadro 
dos SRP conhecidos, ficando agora resolvidos. 


c) as potencialidades práticas são grandes, 
especialmente para os problemas de ordenação 
e selecção de alternativas já que, recorrendo à 
Teoria dos Grafos, podem construir--se funções 
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de potencial que apoiam o decisor com tais 
objectivos. 


À concretização computacional da Hiper - Ordem 
começa por procurar concluir sobre as comparações 
relativamente às quais o decisor não se pronunciou, 
tirando-se partido das comparações conhecidas e das 
condições de cadeia. 


Seguidamente, admitem-se níveis constantes para 
(qg,p) no Qasem análise de detectar eventuais 
heterogeneidades de sensibilidade dó decisor o que 
não é de menor interesse prático. Posteriormente 
exploram-se as ordenações e selecções possiveis de 
estabelecer numa base uni- e multi-critério 
desenvolvendo-se um processo interactivo com o 
decisor 


FIGURA 5 - Exemplo de comparações que exigem q 
variável. 
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Das aplicações feitas no domínio da selecção de 
propostas em concursos Públicos de Empreendimen- 
tos de Engenharia, recolhe-se um exemplo do grafo de 
comparações obtido e que se inclui na Fig.6. 


——e Relações de comparação estabelecidas pelo decisor. 


-— — Relações de comparação não estabelecidas pelo 


decisor mas que são determináveis pelas condições 


de cadela da MHiper-Ordem como sendo P 


FIGURA 6 - Exemplo do grafo de comparações 


5. CONSIDERAÇÕES FINAIS 


Apresentaram-se neste artigo linhas de 
desenvolvimento teórico e metodológico que têm 
contribuido para a racionalização dos problemas de 
decisão. Estas contribuições são uma ajuda valiosa 
para muitos dos problemas complexos vividos pela 
Engenharia moderna o que não significa que a sua 
aceitação não conheça resistência e desconfianças. 
Resistências suscitadas pela inovação que lhes é 
inerente, desconfianças resultantes de se penetrar no 
cerne dos processos de escolha sempre resguardados e 
protegidos ciosamente por cada decisor. E, todavia, 
para além dos diversos aspectos técnicos já referidos, 
julga-se que os métodos e modelos apresentados já 
vêm desempenhando papeis progressivamente 
importantes, em especial no que concerne às funções 
seguintes : 


a) Função semântica explicitando a escassez 
ou O excesso, O rigor ou a imprecisão das 
representações (modelos) do real estudado; 


b) Função sintáctica, estabelecendo sistemas e 
quadros lógicos para o tratamento dos 
problemas representados; 


c) Função integradora de contribuições 
especializadas disciplinares distintas tendo 
em vista a resolução dos problemas em estudo; 


d) Função axiológica, revelando os sistemas de 
valores ou critérios de escolha formulados; 


e) Função interpelativa. contribuindo, através 
de modos interactivos, para a clarificação dos 
espaços menos racionalizados do próprio 
decisor; 


f) Função inventiva, alargando os espaços das 
alternativas possiveis através das 
decomposições e reconstruções morfológicas 
mais adequadas; 


g) Função harmonizadora, garantindo a 
estabilidade de principios e de abordagens nas 
diversas componentes do processo de decisão 
ou detectando incoêrencias e desvios aos 
quadros racionais escolhidos; 


h) Função ética, explicitando a tensão própria 
da relação entre acção empreendida e decisão 
escolhida confirmando-se o significado de 
racionalização desenvolvida ou, pelo 
contrário, revelando-se a subversão parcial ou 
completa do sistema de valores ou critério 
formulados, . 
Espera-se assim, através deste artigo, haver 
contribuido para o esclarecimento dum domínio de 
desenvolvimento cientifico que, segundo tudo indica, 
influenciará, cada vez mais, a forma de pensar e 
realizar a moderna Engenharia. 
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OS CRISTAIS LÍQUIDOS. 
SUAS PROPRIEDADES E APLICAÇÕES 
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Resumo 


Neste trabalho são descritas algumas das fases 
líquido cristalinas conhecidas e são apresentados os 
métodos usualmente utilizados para a identificação e 
caracterização das mesmas. Igualmente se 
apresentam algumas das aplicações, com interesse 
tecnológico, dos cristais líquidos. 


1. INTRODUÇÃO 


Algumas substâncias ao passarem do estado sólido 
ao estado líquido apresentam uma ou várias fases 
intermédias (designadas mesofases) com 
propriedades mecânicas idênticas às dos líquidos mas 
com propriedades ópticas anisótropas próprias dos 
cristais. Assim, as fases líquidas intermédias destas 
substâncias quando observadas num microscópio 
entre polarizadores cruzados, revelam magníficas 
texturas características da anisotropia do meio, 
extinguindo-se a imagem se a substância transitar 
para a fase líquida isótropa. Às substâncias que 
exibem estas fases intermédias foi dado o nome de 
cristais líquidos. 

A obtenção de mesofases pode conseguir-se por 
variação da temperatura (mesomorfismo 
termotrópico), ou por influência de um solvente 
(mesomorfismo liotrópico). 

O primeiro liotrópico conhecido foi obsevado 
casualmente pelo biólogo Virchow |1| em 1854 quan- 


Lisboa 
* - Secção de Ciência dos Materiais(FCT-UNL) Monte da Caparica 


do estudava a estrutura de soluções de mielina em 
água usando um microscópio. O primeiro 
termotrópico foi descoberto em 1888 pelo botânico 
Reinitzer |2] que ao sintetizar o benzoato de 
colesterilo verificou que esta substância apresentava 
"dois pontos de fusão” com características distintas; 
a 145,5C a substância fundia e apresentava um 
aspecto turvo mas a 178,5ºC tornava-se transparente. 
Entre aquelas temperaturas o benzoato de colesterilo 
exibia uma fase colestérica liquída cristalina 
tornando-se liquido isótropo a partir de 178,5ºC. 

Após estas descobertas históricas muito se avançou 
no domínio da síntese química e conhecimento das 
propriedades físicas de novos destes materiais. Nos 
nossos dias são conhecidas múltiplas fases líquidas 
cristalinas, com características estruturais distintas, 
algumas das quais com aplicações imediatas na 
indústria electrónica, na biologia na medicina e 
noutros campos. 


2.CRISTAIS LIQUIDOS LIOTRÓPICOS 


Os cristais líquidos liotrópicos são obtidos a partir 
de determinadas substâncias misturando-as com um 
solvente em determinadas condições de concentração 
pressão e temperatura. Os liotrópicos mais correntes 
são soluções de moléculas anfifílicas em água. Os 
compostos anfifílicos têm uma cabeça polar 
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os solúvel 


Fase Hexagonal Fase Lamelar 


CH, insolúvel 

Pá Fase Micelar 
HL 

" 

FIGURA 1 - 


geralmente solúvel em água e uma cauda orgânica 
insolúvel. 

Para uma dada temperatura e pressão as 
moléculas destas substâncias exibem diferentes 
estados de agregação consoante a concentração do 
solvente. Claro está que as variações de pressão e 
temperatura podem igualmente provocar mudanças 
do estado de agregação. A maior parte dos sabões e 
detergentes dão origem a estruturas liotrópicas 
quando misturados com água. Na figura 1 pode 
ver-se um exemplo de algumas estruturas que é 
possível obter quando se dissolve um sabão em 
águal3|. 


Refira-se que o conhecimento da física dos cristais 
líquidos liotrópicos tem grande importância em 
biologia. Como exemplo citamos o caso das 
membranas celulares que são estruturas formadas 
essencialmente por lípidos, proteinas e água e são 
sistemas liotrópicos. 


3. CRISTAIS LÍQUIDOS 
TERMOTRÓPICOS 


Vamos referirmo-nos sobretudo às estruturas e aos 
estudos de caracterização de mesofases não 
macromoleculares termotrópicas, sendo de notar que 
os materiais que as exibem são obtidos por síntese em 
química orgânica. Tendo em conta a forma das 
moléculas que os constituem os cristais líquidos 
termotrópicos podem ser subdivididos em calamíticos 
e discóticos ( ver quadro |). Os calamíticos são 
constituídos por moléculas alongadas possuindo uma 
parte central relativamente rígida chamada corpo e 
cadeias alifáticas flexíveis nas extremidades 
(exemplo na fig.2). Nos discóticos as moléculas têm a 
forma de discos tendo igualmente um corpo central ao 
qual estão ligadas várias cadeias alifáticas que têm 
tendência a dispor-se no plano do corpo (exemplo na 
fig.3). 


concentração de 


sabão 


Algumas das estruturas liotrópicas que se obtêm quando se mistura um sabão em água.|3]. 


CRISTAIS LÍQUIDOS TERMOTRÓPICOS | 


MACROMOLECULARES | POLMEROS LC) MÃO MACÇADOMOLECULARES 


| CALAMITICOS DISCOTICOS | 


MEMATICOS ESMÉTICOS 


[io] [sado | [tum tos | [Sasita se tata 


COLESERICOS 


[rea pa 


Quadro 1 - Classificação dos cristais liquidos Lermotrópicos 


cut —+O- "= O) es cen AO) cor 


L86"€ L89€ 
Ce + Sa -— -+.- 1 


FIGURA 2 - Molécula calamítica do cristal líquido TB10OA e 


esquema de transições de fase |5]. 


Ra OCOC Ha 


FIGURA 3 


trifenileno e esquema de transições de fase 6). 


- Molécula discótica de um hexa-n-alkanoato de 


Os cristais líquidos calamíticos podem apresentar 
diferentes tipos de polimorfismos exibindo 
nomeadamente fases nemáticas, esméticas e 
colestéricas. 

Nas fases nemáticas as moléculas têm em média 
tendência para ( em domínios macroscópios) se 
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alinharem segundo uma direcção n, designada por 
director |7|, encontrando-se os centros de gravidade 
das mesmas distribuídos de forma aleatória (fig.4). 
Esta orientação das moléculas estende-se a longa 
distância podendo obter-se "monocristais" nemáticos 
usando agentes externos como campos eléctricos ou 
campos magnéticos. Us nemáticos são usualmente 
meios uniaxiais* positivos sendo as fases deste tipo 
pouco viscosas (viscosidades típicas da ordem de 
0.1P). Numa escala decrescente de temperaturas a 
fase nemática ocorre habitualmente imediatamente 
a seguir à fase isótropa, embora existam excepções 
devido ao fenómeno da reentrancia descoberto em 
1975 |8). 


FIGURA 4 -Representação esquemática da fase nemática dos 
calamíticos. 


As fases colestéricas são fases nemáticas torcidas 
em que no seu estado normal as moléculas 
dispõem-se segundo configurações helicoidais (fig.5). 
Estas fases podem obter-se com moléculas quirais ou 
ainda por dissolução de moléculas quirais em 
nemáticos. O passo da hélice colestérica é uma 
função de pressão, da temperatura e também de 
campos exteriores nomeadamente campos eléctricos 
e magnéticos |7,9|. São como os nemáticos pouco 


FIGURA 5 - Representação esquemática da fases colestérica dos 
calamíticos 


Co Ras a 


* Recentemente foram sintetizados nemáticos biaxiais 


viscosos mas são uniaxiais negativos, coincidindo o 
eixo óptico com o eixo da hélice. Estas fases são 
ópticamente activas exibindo poderes rotatórios 
milhares de vezes superiores aos da sacarose |10). 


As fases esméticas são muito mais viscosas do que 
as fases nemáticas e têm como características 
fundamental comum o facto de as moléculas se 
disporem em camadas. Consoante o modo de 
orientação e disposição das moléculas dentro das ca- 
madas, a correlação entre as posições das moléculas 
dentro das camadas e de camada para camada e 
ainda a ordem orientacional de ligação |11| entre as 
moléculas, assim se podem considerar diferentes 
tipos de fases esméticas. Tanto quanto é do conheci- 
mento dos autores existem doze tipos principais de 
fases esméticas distintas ( ver quadro 1). Como 
exemplos iremos descrever as fases esmética A e 
esmética B podendo a descrição detalhada dos 
diferentes tipos de fases esméticas ser encontrada na 
referência |12]. 

Os esméticos À ( SA) assemelham-se a nemáticos 
em camadas. Assim as moléculas dispõem-se em 
média perpendicularmente às camadas sendo 
aleatória a distribuição dos centros de gravidade das 
mesmas dentro das camadas ( fig.6). A ordem de 
posição e a ordem orientacional de ligação são a curta 
distância não havendo correlação entre camadas. 


FIGURA 6 - Representação esquemática de fase S,. 


A fase esmética B (Sp) é uma fase quase cristalina. 
As moléculas dispõem-se perpendicularmente às 
camadas em redes hexagonais centradas (fig.7). A 
ordem posicional e a ordem orientacional de ligação 


FIGURA 7 - Representação esquemática de fase Sa. 
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estendem-se a longa distância existindo de um modo 
geral uma forte correlação entre camadas. 


A terminar esta descrição sobre os cristais líquidos 
calamíticos deve referir-se que num mesmo material 
podem observar-se diferentes mesofases líquidas 
cristalinas (ver o esquema de transições de fase 
representado na figura 2). 


Os cristais líquidos de moléculas discóticas são de 
descoberta recente tendo o primeiro destes materiais 
sido sintetizado em 1977 pelo grupo de 
Chandrasekhar |13). Nos nossos dias já são 
conhecidas oito fases distintas deste tipo de materiais 
(ver quadro [). Citaremos como exemplo as fases do 
tipo D, e também os nemáticos. 


i) Nas fases D, as moléculas dispõem-se em colunas 
na perpendicualr aos eixos das mesmas e estas 
formam uma rede hexagonal centrada. Existem dois 
tipos destas fases : a fase Dho/14] em que a distância 
entre moléculas de uma mesma coluna é constante e 
igual para todas as colunas e a fase Dpdl6| em que não 
existe ordem ao longo do eixo das colunas. Na figura 
& está representado como exemplo o esquema 
estrutural da fase Dp,. 


ii) Na fase nemática discótica Np/15] os planos das 
moléculas têm tendência para se colocarem paralelos 
uns aos outros distribuindo-se os centros de 
gravidade das moléculas de forma aleatória(fig.9). 


Tal como os calamíticos, num mesmo composto de 
moléculas discóticas podem detectar-se diferentes 
mesofases, como se evidencia no exemplo de esquema 
de transições de fase representado da figura 3. 


FIGURA 8 - Representação esquemática da fase D,,. 


À terminar este parágrafo convém referir que na 
maior parte das fases líquidas cristalinas as 
moléculas rodam, difundem, oscilam podendo estes 
movimentos ser mais ou menos livres tudo 
dependendo da natureza da fase. Hoje em dia existem 


técnicas experimentais poderosas como a 
Ressonância Magnética Nuclear (R.M.N.), e a 
difracção de neutrões que nos permitem conhecer os 
tipos de movimentos e a escala de tempos em que 
ocorrem para os diferentes tipos de fases|12]. 


FIGURA 9 - Representação esquemática da fase Np 


4. ALGUNS MÉTODOS DE ANÁLISE E 
CARACTERIZAÇÃO DAS FASES 
LÍQUIDAS CRISTALINAS 
MICROCOSPIA POLARIZANTE, 
MISCIBILIDADE, CALORIMETRIA 
DIFERENCIAL E DIFRACÇÃO DOS 
RAIOSX. 


Uma das técnicas mais utilizadas, desde o início do 
século, para a caracterização dos cristais líquidos é a 
da observação de texturas usando um microscópio 
polarizante. Neste campo é de referir um trabalho 
notável de Friedel |16| publicado em 1922 em que são 
estudadas exaustivamente as texturas dos cristais 
líquidos e proposta com base nessa observação uma 
classificação destes materiais. É no entanto 
necessário realçar que a texturas diferentes podem 
ou não corresponder estruturas diferentes (fases 
diferentes) pelo que em caso de dúvida torna-se 
necessário recorrer a outros métodos nomeadamente 
à difracção dos raios X. Contudo a análise em 
microscopia polarizante é extremamente importante, 
sendo o microscópio, acoplado a uma platina de 
aquecimento, um dos instrumentos fundamentais 
dos químicos de sintese orgânica destes materiais. 
Citaremos a seguir algumas das texturas mais 
comummente encontradas em algumas fases líquidas 
cristalinas : 

1) Fase nemática dos calamíticos - as texturas mais 
comuns são a textura mármore e estriada. Na figura 
10 está representado um exemplo de uma textura 
nemática estriada na qual se detecta o aparecimento 
de caudas negras começando em pontos que 
correspondem a singularidades de orientação do 
director na preparação. De cada ponto de 
singularidade partem duas ou quatro caudas negras 
sendo de referir que na preparação e na zona das 
caudas negras as moléculas estão paralelas ou 
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perpendiculares ao polarizador ou analisador que 
estão cruzados. Rodando simultaneamente 
polarizador e analisador as caudas mudam de posição 
mantendo-se os pontos de singularidade. É também 
de referir que as variações de cor correspondem a 
diferentes orientações das moléculas na preparação. 

ii) Fases esméticas - dada a diversidade das fases 
esméticas seria fastidioso enumerar todas as 
texturas [16,17| observáveis neste tipo de fases pelo 
que apresentamos na figura 11 como exemplo uma 
das texturas mais caracteristicas das fases em 
camadas - a textura em leques com cónicas focais. 

ii) Fases de moléculas discóticas - na fase nemática 
dos discóticos a textura mais comummente observa- 


FIGURA 10 - Textura nemática estriada observada ao micróscópio 
óptico com polarizadores cruzados. 


FIGURA 11 - Textura em leques quebrados com cónicas focais 
observada, entre polarizadores cruzados, na fase Sg de um cristal 
líquido, 


da é a estriada |15| tal como na fase nemática dos 
calamíticos. Na fase Dg, as texturas mais observadas 
são germes homeotrópicos e mosaicos [14,18]. Na fase 
Dha são leques quebrados e cónicas focais |13,14]. Na 
figura 12 está representada, como exemplo, uma 
textura em mosaicos observada na fase Dp, de um 
cristal líquido. 


Outro tipo de estudos usados para caracterizar a 
estrutura dos cristais líquidos são os estudos de 
miscibilidade. Estes estudos mostram-se muito 
interessantes na medida em que conhecendo as 
diferentes mesofases de um dado cristal líquido 
poderemos identificar as mesofases de outros. Estes 


FIGURA 12 Textura em mosaicos observada, com microscópio 
polarizante na fase D,, de um cristal líquido. 


estudos baseiam-se na seguinte regra: para cada 
temperatura se dois materiais mesomorfos forem 
miscíveis em todas as proporções, então as mesofases 
são do mesmo tipo. Contudo a recíproca nem sempre 
é verdadeira visto que podem existir substâncias com 
o mesmo tipo de mesofases e estas não serem 
miscíveis. À análise qualitativa imediata da 
miscibilidade de dois compostos pode ser feita usando 
o chamado método de contacto de Kofler [19] que nos 
permite avaliar o comportamento geral da mistura 
utilizando apenas uma única preparação a observar 
ao microscópico. Para isso os dois compostos são 
postos em contacto sobre a lâmina e fundem-se até 
ambos atingirem a fase isótropa. Descendo a 
temperatura pode observar-se a evolução da região 
de contacto na qual por difusão se forma uma zona 
que contém todas as concentrações possíveis que vão 
de uma substância pura à outra também pura. Por 
variação da temperatura detectam-se as transições 
que ocorrem na zona de contacto podendo as mesmas 
ser registadas fotograficamente (exemplo na figura 
13). O método de contacto é obviamente limitado 
uma vez que não se conhece como varia o gradiente 
de concentrações na preparação. 

O estudo e a determinação exacta dos diagramas 
de fase isobáricos em função de concentração e 
temperatura reveste-se de grande importância pois 
as misturas homogéneas dos cristais líquidos não 
apresentam em geral as mesmas propriedades que 
nos compostos puros. Aliás estudos deste tipo têm 
conduzido a algumas importantes descobertas, entre 
as quais citamos o fenómeno da reentrancia já 
referido |8|] e posteriormente detectado em alguns 
compostos puros. Para se obterem diagramas de fase 
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do tipo representado na figura 14 preparam-se 
diferentes misturas homogéneas com concentrações 
conhecida se para cada uma delas detectam-se as 
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FIGURA 13 - Método de contacto. Exemplo de texturas 
observáveis na zona de contacto entre dois cristais líquidos 
distintos (CBOOA/12CN)|20] 
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FIGURA 14 - Diagrama de fase associado a misturas em diferentes 
concentrações dos cristais líquidos a que se refere a figura 13 
(CBOOA/12CN)|20]. 


temperaturas de transição usando o microscópio 
polarizante e/ou a análise calorimétrica diferencial 
(DSC). A análise calorimétrica diferencial é um dos 
métodos mais precisos utilizados para detectar 
transições de fase tanto em compostos puros como em 
misturas. Com esta técnica é, em geral, possível 
determinar as temperaturas de transição e para as 
transições de fase de 12 ordem as entalpias de 
transição (fig.15). 


À técnica experimental por excelência usada para 
determinar a estrura das fases líquidas cristalinas é 
contudo a difracção dos raios X (exemplo na figura 
16). Esta técnica conjugada comas já descritas 
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FIGURA 15 - 
diferencial(DSC) obtido em temperatura crescente. Os dois picos 
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detectados correspondem a uma Lransição cristal esmético B e a 
outra esmético B isótropo |21] 
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FIGURA 16 - Figuras de difracção de Raios X obtidas na fase 
cristalina (1) e na fase esmética A (Il) de um cristal liquido |20| e 
correspondentes microdensitogramas. 


permite na maior parte dos casos definir com clareza 
o tipo de fases que um dado composto exibe. 


5. ALGUMAS APLICAÇÕES DOS 
CRISTAIS LÍQUIDOS 


Seria fastidioso apresentar uma lista de todas as 
aplicações que podem ser dadas aos cristais líquidos. 
Apresentamos como exemplos três. 


5.1- Os colestéricos como detectores 
térmicos 
Como foi referido em 3. os colestéricos têm 


importantes propriedades ópticas. Em especial é de 
referir que quando um feixe luminoso em certos 
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comprimentos de onda, incide sobre um 
“monocristal” colestérico paralelamente ao eixo da 
hélice do material o feixe é decomposto de um modo 
geral em duas componemtes com polarização circular 
(uma esquerda e outra direita) uma das quais é quase 
totalmente reflectida sendo a outra transmitida tudo 
dependendo da concordância ou discordância entre o 
"sentido" da hélice descrita pelo vector luminoso [5 
num dado instante t e o “sentido” natural de torção 
da hélice do colestérico. Para a incidência nas 
condições definidas verifica-se que a reflexão 
máxima ocorre para um comprimento de onda 
Am = p.n sendo p o passo da hélice e o n o indice de 
refracção médio de uma camada de colestérico. De 
um modo geral os colestéricos exibem um aspecto 
colorido quando iluminados com luz branca sendo a 
cor observada função dos comprimentos de onda 
reflectidos (reflexão de Bragg). Como o passo da 
hélice é uma função da temperatura, variando esta 
podem obter-se variações de cor correspondente à 
zona de comprimentos de onda de reflexão máxima. 


FIGURA 17 - Imagem revelada por um colestérico, colocado entre 


um plástico negro e outro transparente, após contacto de uma 
mão|23|. 


Hoje em dia conseguem-se obter colestéricos ou 
misturas de colestéricos |10,22|] que exibem 
mudanças de cor do vermelho ao violeta em domínios 
desejados de temperatura (grandes ou pequenos 
intervalos). Esta propriedade é usada, entre outras 
aplicações, no teste de febre, na detecção do cancro da 
mama e na detecção de sobreaquecimentos em 
máquinas. Como exemplo na figura 17 encontra-se 
representado um depósito de colestérico sobre um 
fundo negro sendo a imagem obtida após contacto 
com a mão. Às variações de cor correspondem a 
variações locais de temperatura. 


5.2- "Displays" alfanuméricos de 
cristais líquidos. TV a cores com 
cristais líquidos. 


Existem várias técnicas para construir os displays 
alfanuméricos a cristal líquido que são usados em 
relógios, calculadoras, letreiros informativos, painéis 
de bordo de automóveis etc. 

Na técnica TN (nemático torcido) que é uma das 
mais conhecidas |24|] usa-se um cristal líquido 
nemático que é introduzido numa célula limitada nos 
topos por dois vidros transparentes cada um deles 
com um depósito condutor também transparente. À 
distância entre os vidros situa-se habitualmente 
entre 10 e 20 micron. Ás faces das placas de vidro 
que ficam em contacto com o cristal líquido são 
previamente preparadas de forma a forçarem que 
junto às mesmas as moléculas do nemático se 
orientem paralelamente aos vidros mas com o 
director rodado de 900 quando se passa de uma placa 
à outra (fig. 18). Deste modo as duas condições 
fronteiras induzem as moléculas do cristal líquido a 
disporem-se, no interior da célula, segundo uma 
configuração helicoidal (fig. 18) semelhante à 
observável expontâneamente nos colestéricos. À 
célula é limitada por dois polarizadores cruzados de 
tal forma que radiação luminosa despolarizada que 
incide sobre um deles abandona-o com o vector 
luminoso E a vibrar segundo a direcção de 
alinhamento das moléculas junto da face de vidro 
mais próximo. Esta direcção de vibração é no entanto 
rodada, devido a estrutura helicoidal adoptada pelo 
cristal líquido, de tal forma que a luz emergente 
encontra o segundo polarizador em condições de 
transmissão total (fig. 18), No nemático utilizado na 
figura 18 a anisotropia dieléctrica é positiva (À £>0) 
e deste modo quando se aplica um campo eléctrico 
entre as placas da célula as moléculas têm tendência 
a alinhar-se segundo o campo (configuração 
homeotrópica). Nestas condições luz linearmente 
polarizada proveniente do 1º polarizador atravessa a 
célula sem que se verifique rotação da direcção de 
polarização da mesma provocando assim o segundo 
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polarizador a extinção (fig. 18). Nestas condições o 
display apresenta-se escuro. 
A descrição apresentada permite com preender o 
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FIGURA 18 - Princípio de funcionamento dos displays de nemático 
torcido (EN). 


princípio físico de funcionamento das células TN 
embora seja necessário esclarecer que na realidade só 
uma das placas de vidro tem uma face inteiramente 
condutora tendo a outra segmentos condutores sobre 
uma base não condutora. Assim, quando se aplica 
tensão entre os segmentos e a base condutora são 
visualizados caracteres escuros sobre um fundo claro 
sendo o contraste assim obtido notável (ver fig.19). 
Este principio é usado em letreiros iluminados 
posteriormente ou não. No caso dos grandes letreiros 
informativos existe iluminação posterior enquanto 
que nos displays dos relógios e outros não existe este 
tipo de iluminação. Nestes na face posterior da célula 
sob o polarizador inferior existe uma superfície que 
reflecte a luz ambiente que a atinge podendo assim 
observar-se os efeitos descritos, 
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FIGURA 19 - Representação esquemática de uma célula real 


numérica de nemático torcido CTN. 


Nos últimos anos começaram a ser introduzidos no 
mercado televisores a cores com ecran a cristal 
líquido de pequena e média dimensão |24). Têm em 
relação aos televisores convencionais de raios 
catódicos a vantagem imediata de ocuparem um 
muito menor volume pois não é necessária a 
espessura considerável absolutamente indispensável 
nos televisores convencionais para aceleração e 
focalização dos electrões. Um ecran a cores de cristal 
liquido é dividido em unidades elementares. Um 
ecran de 20xl5cm tem aproximadamente 100.000 
destas unidades elementares. Cada unidade 
elementar é formada por três células de nemático 
torcido (TN). Uma célula possui um filtro azul, outra 
um filtro verde e a outra um vermelho. Quando 
nenhuma voltagem é aplicada cada célula reflecte a 
cor do filtro que possui; quando é aplicada tensão a 
célula fica preta. As diferentes cores obtêm-se por 
diferentes combinações on-off de cada unidade 
elementar. Na figura 20 está representada, como 
exemplo, uma imagem obtida num televisor a cores 
portátil que usa a técnica TN. 

Tudo aponta para que nos próximos anos a 
tecnologia de multiplexagem e comando electrónico 
se desenvolva de forma a ser possível construir e 
comercializar ecrans deste tipo de grandes 
dimensões. 


FIGURA 20- Imagem obtida num televisor portátil a cristais 
liguidos que usa a técnica TN |24 e |. 


d.3- Nova Tecnologia de Displays a 
Cristal Líquido 


Uma nova tecnologia de displays que faz uso de 
cristais líquidos começou recentemente a ser 
desenvolvida em vários laboratórios nos EUA e 
Japão |25]. Nos displays fabricados usando esta 
tecnologia, a modulação da luz é obtida controlando a 
difusão luminosa causada por bolhas de cristal 
líquido na fase nemática dispersas num polimero 
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transparente. O princípio de funcionamento destes 
displays é baseado na birrefringência da fase 
nemática e no facto de que o campo do director no 
interior das bolhas para certos cristais líquidos e 
certos polímeros apresentar uma estrutura bipolar 
com dois defeitos pontuais em pontos opostos da 
bolha. 

Um display construído segundo esta tecnologia está 
exemplificado na figura 21. As placas de vidro 
contêm o depósito metálico transparente já referido 
em 5.2 que permite aplicar um campo eléctrico 
perpendicular às mesmas. Entre as placas 
encontra-se um polímero transparente com um índice 
de refracção igual ao índice de refracção ordinário do 
nemático contido nas bolhas. Na ausência de 
qualquer campo exterior as estruturas do director no 
interior das diferentes bolhas estão orientadas ao 
acaso, sendo a orientação imposta pela morfologia 
local. Um raio luminoso incidindo normalmente no 
display encontra no seu percurso variações do indice 
de refracção ao passar em bolhas que não estejam 
alinhadas com a normal ao display e é difundido.Para 
bolhas com um diâmetro da ordem do micron o 
display é opaco e apresenta uma cor branca. Quando 
um campo eléctrico suficientemente intenso é 
introduzido no display, por aplicação de voltagem 
entre os eléctrodos, as estruturas do director nas 
diferentes bolhas alinham com o campo eléctrico e os 
raios luminosos passam pelo display sem 
encontrarem variações do índice de refracção e não 


são difundidos. O display torna-se transparente. À 
aplicação do campo eléctrico muda o display de um 
estado de difusão para um estado transparencia. 
Existem basicamente quatro processos de produzir o 
polímero transparente com as bolhas de cristal 
líquido dispersas no seu interior : 
a) Encapsular o nemático por acção de um agente 
quimico dentro de pequenas bolas que são depois 
misturadas com o polímero transparente |26]. 
b) Separação de fases entre o cristal líquido e o 
polímero induzida por polimerização |27). As resinas 
epoxidícas podem ser usadas nesta aplicação.O cris- 
tal líquido é misturado uniformemente com o 
componente À da resina e é feita seguidamente a 
mistura com a componente B obtendo-se 
inicialmente uma mistura homogénea. À medida que 
a polimerização vai ocorrendo o cristal líquido 
torna-se insolúvel e concentra-se em bolhas que 
surgem dispersas no polímero. O processo de 
formação das bolhas pode ser regido quer por 
nucleação e crescimento |28], quer por decomposição 
spinoidal [29]. Controlando a temperatura a que se dá 
a polimerização é possível controlar o diâmetro das 
bolhas assim como a dispersão neste valor. 
c) Separação de fases induzida pela temperatura |27|. 
Usa-se um polímero termoplástico que se torna 
solúvel com o cristal líquido a alta temperatura. 
Quando a temperatura é reduzida o cristal líquido 
torna-se insóluvel e formam-se bolhas no polímero. 
Controlando a taxa de arrefecimento é possível 
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FIGURA 21 - Display de bolhas de cristal líquido dispersas num polímero transparente na situação de difusão (sem campo 


aplicado). 
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controlar o diâmetro das bolhas assim como a 
dispersão neste valor. 

d) Separação de fases induzidas por evaporação |27]. 
Prepara-se uma mistura homogénea de cristal 
líquido e polímero transparente usando um solvente 
apropriado. Este solvente é depois evaporado e o 
cristal líquido torna-se insolúvel concentrando-se em 
bolhas que aparecem uniformemente dispersas no 
polímero. 


Um ponto de particular interesse no 
desenvolvimento destes displays são os tempo on e off 
assim como a voltagem necessária para mudar o 
estado do display de difusão para transparente. 
Valores típicos para os tempos on e off são da ordem 
de alguns milisegundos. O campo eléctrico para 
actuar o display é presentamente da ordem de 
1,5Volt/um |30]. Estes pontos são presentemente foco 
de intensa pesquisa. O contraste nestes displays 
pode ser controlado quer variando a espessura, quer 
ajustando o índice de refracção do polímero de modo a 
ser tão próximo quanto possível do índice ordinário 
do cristal líquido. Na figura 22 podem ver-se dois 
displays reais e idênticos, colocados lado a lado, 
fabricados segundo o processo b). O da direita 
encontra-se no estade de transparência por aplicação 
de uma tensão de valor eficaz igual a 15V, sendo 
assim possível visualizar o relógio colocado 
posteriormento, e o da esquerda no estado de 
dispersão sem tensão aplicada. 


Estes displays apresentam um campo bastante 
vasto de possíveis aplicações, que não é acessível aos 
displays de nemático torcido (TN). A técnica de 


FIGURA 22 - Displays de bolhas de cristal líquido dispersas num 
polímero. O display da direita encontra-se no estado de 


transparência e o da esquerda encontra-se no estado de dispersão. 


fabrico é bastante simples e presta-se facilmente à 
produção de displays de grande área. Uma das 
aplicações é em sinalização de baixo conteúdo 
informativo : sinais de tráfego, anúncios para 
estabelecimentos, letreiros informativos, etc.. Estes 
materiais têm também uma grande potencialidade 
na construção de edifícios, pois é possível obter 
vidros que mudam de opacos para transparentes por 
simples aplicação de uma pequena voltagem. Um 
outro campo de aplicação com grandes possibilidades 
é a indústria de automóveis, quer para as janelas, 
quer para os mostradores. 

Em aplicações de alto conteúdo informativo, estes 
materiais requerem uma matriz activa pois possuem 
um campo eléctrico crítico mal defenido. Usando uma 
matriz activa estes displays podem ser usados em 
sistemas de projecção de imagem e também como 
écrans planos para televisão. 
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A SÍNTESE DUM PROCESSO 
QUÍMICO INDUSTRIAL, 
UTILIZANDO A TECNOLOGIA DO 
PONTO DE ESTRANGULAMENTO 


CLEMENTE PEDRO NUNES(1) E BODO LINNHOFF (2) 


SUMÁRIO 


As técnicas de análise e de síntese destinadas ao 
desenvolvimento e ao projecto dum processo químico 
industrial evoluiram significativamente nos últimos 
trinta a quarenta anos. Us processos industriais têm 
de ser hoje mais sóbrios em termos de capital e, mais 
flexíveis e eficientes em termos de energia e 
rendimento. 


Os avanços nos conceitos teóricos relativos ao 
projecto químico resultaram não só numa melhor 
simulação e modelação como também, numa melhor 
compreensão de princípios que determinam a Síntese 
do Processo. Em particular a Tecnologia do Ponto de 
Estrangulamento (Pinch Technology) conduziu a 
importantes poupanças ( uma “alteração 
qualitativa") na prática industrial através duma 
mais correcta integração do processo. 


O presente trabalho analisa a evolução dos conceitos 
de síntese no contexto do Ensino de Engenharia 
Química até meados dos anos 60, bem como os 
desenvolvimentos posteriores com particular enfase 
na síntese do processo e na tecnologia do Ponto de 
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Estrangulamento. É referida a experiência verificada 
nalgumas importantes Escolas de Engenharia 
Química e, em particular, a alteração recentemente 
verificada no Instituto Superior Técnico, em Lisboa. 


O ensino dos novos conceitos de síntese é ilustrado 
com um exemplo simples dum diagrama de fluxos 
dum processo ao qual vão sendo aplicados de forma 
concreta esses conceitos culminando com, a 
tecnologia do Ponto de Estrangulamento; assim 
sepodem obter diagramas de fluxos alternativos para 
o mesmo processo que, conduzem a notáveis reduções 
de consumo de utilidades. 


Finalmente o presente artigo aborda a dificuldade 
densinar "Projecto" (o "Projecto" é criativo e 
imaginativo e, como tal só parcialmente pode ser 
ensinado). Sugere-se que o tema do desenvolvimento 
global do processo oferece uma excelente 
oportunidade para ensinar Projecto em engenharia 
química e, até certo ponto, permite às Universidades, 
formar engenheiros com alguma experiência. 


a. Departamento de Engenharia Quimica Instituto Superior Técnico,Lisboa, Portugal. 


a Department of Chemical Engineering University of Manchester . Institute of Science and Technology, Manchester, Englnico. 
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SUMMARY 


Analysis and synthesis techiques for the 
development and design of the overall chemical 
process have undergone significant changes in the 
last thirty to forty years. Processes today have to be 
leaner in terms of capital, more flexible, and more 
efficient in terms of energy and yield. 


Advances in theory have resulted in better modeling 
and simulation as well as improved understandig of 
the principles governingProcess Synthesis. In 
particular Pinch Technology ( Tecnologia do Ponto de 
Estrangulamento) has resulted in significanth 
pratical savings (a "Step change") through better 
process integration, 


This paper discusses the evolution of synthesis 
concepts within Chemical Engineering Teaching to 
the mid 60's and the developments since with 
particular emphasis on process synthesis and Pinch 
technology. It is described the experience of a 
number of important Chemical Engineering Schools 
and, specially the change recently registered in 
Instituto Superior Técnico, Lisbon. 


The teaching of the new synthesis concepts is 
ilustrated through a simple example of a process flow 
diagram to which it is successively apllied those 
concepts, up to the Pinch technology which lead to 
remarquable improvements in the utility 
consumption, 


Finally, the paper discusses the difficulty of teaching 
“ design” (design is creative and imaginative and, as 
consequence it can only be taught in a limited way). 


It is suggested that the subject of overall process 
development offers a splendid chance to teach 
Chemical Engineering design and, up to a point, 
offers a chance to universities Lo Lurn out 
experienced engineers. 


1. OS CONCEITOS DE "ANÁLISE" E 
DE "SINTESE" NO DESENVOLVI- 
MENTO DUM PROCESSO DE 
ENGENHARIA QUÍMICA. 


O estudo do desenvolvimento dum processo químico 
industrial requer conhecimentos e tecnologias de 
síntese e de análise, necessários à elaboração do 
respectivo projecto. 


E o projecto químico é, como referiu Asimov, o 
principal objectivo e o aspecto mais marcante da 
actividade do engenheiro químico. 


1.1 A "Análise" 


A análise é um termo familiar tanto para os 
estudantes como para os diplomados de engenharia 
química; geralmente a análise processa-se pela 
decomposição do sistema nas sua partes ( ou 
equipamentos) componentes, estudando-se em 
profundidade cada um desses equipamentos em 
separado. 


Desde a 28 guerra mundial o ensino da engenharia 
química acentuou claramente a sua componente de 
análise dedicando particular atenção ao estudo das 
operações unitárias, fenómenos de transporte, 
instrumentação e controlo de processos e engenharia 
das reacções químicas. 


1.2 A "Sintese” 


Pelo contrário não houve, até muito recentemente, 
um desenvolvimento paralelo das disciplinas que 
estudam a sintese do processo, que se pode definir 
como "a combinação e articulação inteligente dos 
diversos conceitos e conhecimentos fornecidos pela 
análise, num conjunto coerente que irá promover a 
realização dum determinado objectivo industrial, de 
forma técnica e economicamente optimizada”. 


No ciclo interactivo que leva à realização do projecto 
químico a síntese precede a análise mas a análise 
sugere normalmente mudanças à sintese. 


Assim, durante o desenvolvimento do projecto, a sua 
optimização depende tanto da síntese como da 
análise. 


2. A IMPORTÂNCIA DADA AO 
ESTUDO DA SÍNTESE DAS VÁRIAS 
ETAPAS DUM PROCESSO, NOS 
CURRICULA DO ENSINO DA 
ENGENHARIA QUÍMICA. 


O ensino da engenharia química, nomeadamente em 
Portugal, tem registado um apreciável 
desenvolvimento no estudo da análise do processo, 
com o tratamento quantificado de vários 
equipamentos e fenómenos, enquanto que o ensino da 
síntese tem permanecido até agora restringido a 
técnicas indutivas relativamente pouco 
desenvolvidas. 
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É de salientar que a síntese, ao nível do projecto 
químico, necessita de considerável soma de 
conhecimentos e de experiência dos princípios em 
que se baseia a análise, antes de poder ser aplicada 
com eficiência. 


Neste contexto, o estudo da síntese do processo 
químico havia merecido até à ultima reestruturação 
curricular, uma atenção relativamente pouco 
desenvolvida no Currículo do Curso de Engenharia 
Química do IST; o estudo das técnicas de síntese 
centrava-se na elaboração de balanços globais de 
massas e de energia dum processo, do fluxo de 
informação das variáveis do processo no flowsheet e, 
na optimização de problemas simplificados. 


Internacionalmente e, excepto para o tratamento 
isolado da síntese de redes de permutadores de calor 
e de sequência de colunas de destilação para 
misturas multicomponentes, a preocupação com o 
desenvolvimento de técnicas formais para a síntese 
do projecto químico apenas teve início a partir de 
1968, nos Estados Unidas da América. 


Recorde-se que o objectivo das tecnologias de síntese 
em engenharia química, é o de "equacionar e 
seleccionar, duma forma sistemática, as múltiplas 
alternativas para articular entre si as várias etapas 
do processo químico”. 


Assim há que enfrentar questões como : 


a) Que tipos de operações unitárias de 
separação se devem utilizar e que articulação 
deve ser estabelecida entre elas, no 
tratamento dos produtos duma reacção de 
forma a satisfazer as exigências do mercado ? 


b)JQual a sequência de separações que se deve 
utilizar para, a partir de uma mistura 
multicomponente, se obterem os vários 
componentes puros? 


c)JEm que condições é que será vantajoso um 
processo incluir o aproveitamento de 
subprodutos? 


dJQual o número e sequência de permutadores 
de calor que optimiza um processo ? 


eJQuais as quantidades mínimas de entalpia a 
fornecer (aquecimento) e, a retirar 
(refrigeração) do processo, de forma a 
optimizá-lo ? 


A este tipo de questões a prática da engenharia 
química internacional ( o "state of the art" como é 
correntemente designado) respondia até 1968 com 


processos industriais específicos normalmente 
patenteados pelas mais importantes empresas de 
“engineering” e baseadas em meia dúzia de regras 
heurísticas. 


Estas regras heurísticas visavam sobretudo 
simplificar o tratamento meramente matemático das 
muitas alternativas que se apresentam num 
problema de síntese mesmo relativamente simples. 
Registe-se por exemplo que para separar por colunas 
de destilação uma mistura com oito componentes, nos 
seus constituintes, há 429 alternativas. 


Assim, relativamente à sequência das operações de 
destilação a primeira regra que foi apresentada, 
estabelecia que os componentes devem ser removidos 
um a um , sob a forma de destilados, a partir dos mais 
voláteis e, através duma série de colunas de 
destilação. Esta é a regra de sequência directa de 
volatilidades decrescentes que, geralmente conduz a 
uma redução da pressão de trabalho. 


Mas já em 1947, Lockhart “apresentava resultados 
aplicados ao tratamento de gasolina natural, em que 
a sequência directa deixava de conduzir à situação 
óptima quando o componente menos volátil era o 
constituinte principal da alimentação. Neste caso 
sugeria-se uma sequência indirecta. 


Curiosamente em 1957, Herbert”, numa antecipação 
da articulação de critérios de balanços de massa e de 
energia, concluia que a sequência preferencial deverá 
ser aquela que requer uma quantidade mínima de 
entalpia total a fornecer pois, este é o factor 
dominante na avaliação económica. 


Daqui deduziu duas regras : 
- a sequência deve conduzir a que em cada 
coluna o destilado e o resíduo devem ter 
quantidades sensivelmente iguais e; 
- as separações mais difíceis de levar a cabo 
devem ser efectuadas com um mínimo de 
alimentação, ou seja, colocadas o mais tarde 
possível na separação. 
Mais tarde as técnicas de integração de processos 
foram tratados de forma sistemática sob a orientação 
de D. Rudd do Departamento de Engenharia Química 
da Universidade de Winsconsin, tendo vários dos 
seus colaboradores aperfeiçoado estas regras 
heurísticas. 


Assim na sua tese de doutoramento nesta 
Universidade, Gary Powers em 1971”, sugeriu a 
utilização de quatro regras para determinar a 
sequência óptima de colunas de destilação : 


*em princípio a sequência directa deverá ser 
favorecida; 
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*se um componente for predominante numa dada 
alimentação, removê-lo primeiro, 


*deixar a(s) separação(ões) mais dificeis para o final 
da sequência; 


*remover todos os componentes corrosivos, 
eventualmente presentes, no início da sequência. 


3. UM EXEMPLO SIMPLES DUM 
DIAGRAMA DE PROCESSOS. 
DETERMINAÇÃO DOS RESPECTIVOS 
CONSUMOS ENERGÉTICOS. 


Pretende-se estabelecer um diagrama de fluxos que 
melhor corresponda a um processo constituido 
essencialmente por um reactor em que dois reagentes 
AeB reagem entre si para darem o composto C. 


O reagente B é o reagente limitante; mas a sua 
conversão não é total; por razões económicas 
pretende-se recircular Lodo o reagente À não 
reagido. 


Os composto A, B, e C encontram-se no estado líquido 
à temperatura ambiente, tendo a seguinte ordem de 
volatilidades: 

A>B>C 
não havendo indicação de qualquer deles ser 
corrosivo ou que exija qualquer tipo especial de 
material, 


A corrente de saída do reactor constituída por 
A+B+C pode assim ser separada nos seus 
componentes através de duas colunas de destilação 
colocadas em sequência directa, Lanto mais que o 
componente À é o que se encontra em maior 
percentagem dos três. 


A 12 . coluna de destilação separa o componente 
mais volátil, A, como destilado enquanto que a 28 , 
coluna provoca a separação dos dois componentes 
menos voláteis, B como destilado e € como resíduo. 


Sabe-se que as condições cinéticas da reacção levam 
a que a alimentação do reactor se processe a 1700C, 
sendo a reacção fortemente exotérmica, todavia O 
correcto funcionamento da 12 . coluna de separação 
exige uma alimentação a 1800C€, a fim de garantir a 
necessária pureza do destilado A . 


De acordo com estes dados foi elaborado o diagrama 
processual | apresentado a seguir, a pressão é à 
partida um invariante tanto no reactor como nas 
colunas da destilação. 


Neste processo as correntes que têm que ser sujeitas 
a fornecimentos ou a retiradas de entalpia são as 
correntes [1], [2], [3], |4], [5] e [6] devidamente 
assinaladas no Diagrama Processual I; 


FIGURA 1 


Para cada uma desta correntes é calculado, de acordo 
com o Processo, o respectivo caudal mássico por hora 
e determina-se, a partir de dados tabelados, o 
respectivo calor especifico médio no intervalo de 
temperatura em estudo. 

O produto destes dois valores 


Fx(C) 


p médio 


=Cupacidade Calorífica da Corrente 


traduz a capacidade calorífica da corrente e, neste 
exemplo, tem por unidades: 


Kcal/0C.h 


Note-se também que esta capacidade calorífica média 
inclui todos os factores que alteram o conteúdo 
entálpico da corrente devido à mudança de 
temperatura indicada e, incluindo logicamente os 
calores da vaporização, de condensação ou outros. 


HH 


Cap Calda Corrente = 
AT 


Os dados relativos às correntes a que é preciso 
fornecer ou retirar entalpia neste diagrama de 
processos são apresentados no Quadro 1. 


QUADRO 1 


Corrente Ea Cpm 
aii Kal h 


CO ET E 
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Os 6 permutadores de calor que é necessário instalar 
para satisfazer este Diagrama de Processo, são todos 
em contra corrente por serem os do tipo mais 
eficiente. 


Considera-se também para simplificar, que as 
respectivas áreas de permuta entre correntes foram 
calculadas de forma a que o ATmin fosse igual para 
todos eles: 


DE an ane 


| 


- Calculo do consumo de utilidades deste processo 


De acordo com o Diagrama Processual | há três 
permutadores em que o processo recebe entalpia do 
exterior (utilidade quente) e, três permutadores em 
que é retirada entalpia ao processo (utilidade fria). 


Tradicionalmente a forma mais simples e 
aparentemente mais expedita de resolver o problema 
consistia no fornecimento directo de utilidades 
sempre que tal fosse requerido pelo processo, o que 
tinha vantagens no controlo independente e 
autónomo do processo ao nível das operações 
unitárias e, como tal se manteve até à década de 70. 


É este o princípio que conduz ao Diagrama 
Processual I que, retrata o primado da análise sobre a 
sintese . 


*A entalpia requerida da utilidade quente (vapor) 
corresponde pois a: 


(120 x 2,5) + (20x 25) + (10x 16)= 960 x 103Kcal/h 


*A entalpia que a utilidade fria necessite de retirar 
do processo corresponde pois a: 


(210x 2,5) + (20x 20) + (14x 10) = 1065 x 103Kcal/h 


3.1 Optimização heurística do consumo de 
utilidades deste processo. 


A primeira abordagem para a transformação dum 
diagrama de processos, com vista a reduzir o seu 
consumo de utilidades, baseava-se numa análise 
“local” das correntes quentes e frias dum mesmo 
equipamento ou, de equipamentos adjacentes; assim 
procurava-se a suboptimização de zonas de processo 
que, posteriormente, iriam ser integradas entre si, 
sempre que tal fosse possível. 


Neste diagrama processual a suboptimização local 
mais óbvia envolve o pré-aquecimento da corrente [1] 
através da entalpia excedentária somente na 
correntel2). 


Desta forma no permutador Pl o consumo de vapor 
pode ser eliminado por substituição pela entalpia a 
fornecer pela corrente 2, em que esta corrente irá 
abandonar Pl à temperatura de : 


120x2,5=(390-T)x 2,5 => T=2700C 


Assim esta corrente [2] terá ainda de ser arrefecida 
com utilidade fria no permutador P2, a fim de atingir 
os 180º0C exigidos pelo processo à entrada da 18 
coluna de destilação 


Obtém-se assim uma redução significativa do 
consumo de utilidades: 
*entalpia requerida da utilidade quente(vapor): 


0+(20x 25) +U0x 16)=660x 108Kcal/h 


*entalpia que a utilidade fria (água de refrigeração) 
necessita de retirar ao processo: 

(90 x 2,5) + (20x 20) + (10X 14) = 765 x 103Kcal/h 
Houve assim uma redução de 300x 103Kcal/h tanto 
no consumo da utilidade quente como no consumo da 
utilidade fria. É de salientar que em consequência do 
1º Princípio da Termodinâmica, uma redução no 
consumo da utilidade quente terá forçosamente, 
como bónus, uma redução equivalente no consumo da 
utilidade fria. 


Dever-se-ã agora verificar se exisle alguma 
possibilidade adicional de poupança de vapor; 
recordamos todavia que as poupanças a efectuar 
deverão estar submetidas a dois princípios gerais: 


*o (AT) min de todos os permutadores será de 20ºC; 


*o número de permutadores não poderá ser superior 
ao indicado no Diagrama Processual | que deverá ser 
considerado como Diagrama de referência; 
neutralizando-se assim a variável 
suplementar"capital investido”. 


Assim, a corrente[2] a 2700C, que é a única corrente 
quente disponível atendendo ao seu conteúdo 
entálpico e ao (AT)min, apenas pode aquecer a 
corrente [3] mas de forma totalmente insuficiente, 
continuando a ser necessário um aquecimento 
suplementar com vapor à corrente [3] e um 
arrefecimento suplementar com água de refrigeração 
à corrente [2]. 


Conclui-se assim que, com o mesmo número de 
permutadores apresentado no Diagrama de 
referência, se encontra esgotada aparentemente a 
possibilidade de se reduzir o consumo de utilidades 
neste processo. 

O Diagrama Processual II representa esta nova 
situação, localmente optimizada. 
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4. AS TECNICAS DE OPTIMIZAÇÃO 
DA SINTESE DE PROCESSOS 
DESENVOLVIDOS POR RUDD, MAS- 
-SO, POWERS, SIIROLA E ICHIKAWA. 


Havendo a consciência nalgumas escolas de 
engenharia química, mais voltadas para a síntese do 
Processo, de que a integração de processos se 
encontrava ainda numa fase casuística procurou-se, 
a partir de meados da década de 60, obter diagramas 
processuais (flowsheets)que optimizassem os 
consumos das utilidades energéticas. 


Com esse objectivo, sucessivos investigadores 
preocuparam-se sobretudo com a forma de proceder 
às  interligações entre  permutadores 
(Heat-Exchangers Network - HEN). 


Apenas a partir de 1968, sobretudo por parte do 
grupo do | + D de Dale Rudd, se começaram a deduzir 

regras sistemáticas para simplificar as combinações 

possíveis de permutas de calor entre correntes 


FIGURA 2 
quentes e frias num processo; de facto como Rudd e 


Masso observaram em 1969, um exemplo 
extremamente simples contendo apenas duas 
correntes quentes e duas frias, pode conduzir a 4 200 
combinações possiveis de processar em conjunto estas 
4 correntes. 


Rudd e Masso '!'consideraram nomeadamente as 
seguintes regras : 
*escolher combinações que resultem numa 
recuperação máxima de entalpia disponível numa 
corrente quente; 


*recorrer apenas às utilidades exteriores quando 
estiverem esgotadas as possibilidades de recolher 
entalpia a partir do próprio processo. 


Mais tarde os mesmo autores e K, Lee!t%), aplicaram a 
estratégia de optimização conhecida por 
programação dinâmica, através da qual melhoraram 
significativamente as soluções alcançadas desde que 
os diagramas das redes de permutadores fossem 
totalmente acíclicos. 


Nesta fase a contribuição do Dep. de Engenharia 
Industrial do Instituto de Tecnologia de Tóquio em 
colaboração com a Chiyoda Chemical Engineering 
Co., foi de grande importância. 


Poderosos métodos matemáticos de optimização 
apresentados por Ichikawa, Kobayashi e outros , em 
1971 (6) permitiram estabelecer regras que se 
revelaram capazes de melhorar significativamente 
alguns processos comercialmente disponíveis para a 
indústria: 

- Às trocas de calor entre correntes do processo devem 
ser efectuadas em permutadores de calor em 
contracorrente; 

- Às correntes quentes (e frias) devem ser 
combinadas em ordem decrescente ( e crescente) e 
consecutiva das respectivas temperaturas, 

- À temperatura de saída dum permutador, duma 
corrente quente que está a ser arrefecida, não deverá 
ser nunca inferior à temperatura de saida da 
corrente fria sequencialmente mais quente. 


Estas regras foram em grande parte aceites no livro 
“Process Synthesis" de Rudd, Powers e Siirola (7) que, 
todavia chamou a atenção para o problema destas 
regras poderem não ser consistentes com os (AT) min 
dos diversos permutadores. 


Saliente-se ainda que as regras deduzidas por 
Ichikawa e outros, se baseavam em exemplos onde 
tinham sido introduzidas simplificações que, no seu 
conjunto, são bastante drásticas: 


*O custo dos sistemas de permutadores de calor ser 
uma função linear das áreas de transferência; 

*Os permutadores serem de camisa tubular e, em 
contracorrente; 

*Considerar-se não haver mudanças de estado nas 
correntes do processo; 

*Admitirem-se valores iguais para os coeficientes 
globais de transferência de calor; em todos os 
permutadores. 


Apesar dos seus inegáveis méritos, eslas regras 
aperfeiçoadas por Ichikawa, não garantem ainda, 
como referiram mais tarde Rudd e outros '8), que se 
atinja uma situação óptima g | | processo 
pois, muito simplesmente, não se podia ainda definir 
o valor óptimo a atingir, em termos de consumo de 
utilidades, para o processo global. 


5. UMA NOVA ESTRATÉGIA PARA A 
OBTENÇÃO DE DIAGRAMAS DE 
PROCESSOS A QUE CORRESPONDA 
UMA OPTIMIZAÇÃO DO CONSUMO 
DE UTILIDADES ENERGÉTICAS. 


Vimos nos pontos anteriores quais os esforços 
desenvolvidos até 1975, para se poderem elaborar 
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diagramas de fluxos dos processos globais, de forma a 
optimizar os consumos das utilidades energéticas. 


Todavia, como acabámos de ver, desconhecia-se 
ainda qual o Consumo Mínimo Teórico de Utilidades 
Quentes e Frias de um determinado processo; como 
se continuava ainda a desconhecer o efeito 
correlacionado entre consumos de utilidades e 
investimentos em permutadores de calor. 


E, noutra perspectiva, não tinha ainda sido 
apresentada uma ligação sólida e cientificamente 
fundamentada entre os princípios fundamentais da 
termodinâmica e, o design dos vários diagramas de 
fluxos que se vão obtendo com vista à optimização do 
processo. 


- Como definir , de forma precisa e com base 
científica, a quantidade mínima de utilidades 


quentes e frias e, as respectivas características, que é 
necessário consumir num determinado processo? 


O conceito de Balanço Entálpico do processo, baseado 
no 1º Princípio da Termodinâmica, não resolve o 
problema pois não entra em linha de conta com a 
necessidade simultânea de fornecimento de entalpia 
e de retirada da entalpia do processo, a vários níveis 
diferentes de temperatura. 


Neste ponto e, como se sublinhou no início, a ciência 
da engenharia química avançara muito no aspecto 
analítico (determinação de calores específicos, 
latentes, de misturas, etc.) mas muito pouco no 
aspecto da síntese. 

Valerá a pena referir a este respeito que o livro " 
Introduction to Chemical Engineering" de Edward 
Thompson e William Ceckler do Departamento de 
Engenharia Química da Universidade de Maine e 
publicado em 1977, apenas refere o balanço entálpico 
do processo, sem fazer qualquer referência ao efeito 
do traçado do diagrama de fluxos na determinação do 
consumo energético. 


- Era pois necessário avançar. Mas como? 


d.1 A determinação do ponto de 
estrangulamento ("pinch") duma cascata de 
calor dum processo. Quantidade mínima de 
utilidades teoricamente necessárias. 


De acordo com o 2º Princípio da Termodinâmica o 
calor transfere-se duma fonte quente para uma fonte 
fria. Para um engenheiro de processo a 
disponibilidade de calor só tem significado se, esse 
calor estiver em condições térmicas de ser transferido 
para a corrente que se pretende aquecer. 


Utilizando o conceito de Capacidade Calorífica de 


uma Corrente já atrás definido, podemos agrupar, 


para cada intervalo de temperatura em que só 
participam o mesmo conjunto de correntes, por um 
lado, todas as correntes quentes (ou seja que cedem 
calor) e, por outro lado, todas as correntes frias (ou 
seja as que recebem calor). 

Admitindo como fez Ichikawa t6) e, o livro "Process 
Synthesis" de Rudd e outros ( que, os permutadores 
utilizados são de camisa tubular em contra corrente 
então, em termos teóricos, as Correntes Globais 
Quente e Fria podem permutar calor até ao limite 
imposto pelo (AT) min dos permutadores de calor. 


Para facilitarmos a representação gráfica da Cascata 
de Calor, estabelecemos uma nova referência de 
temperatura, T”, em que: 


*para as temperaturas das correntes quentes: 
T=T1/2(A7) min 

*para as temperaturas das correntes frias : 
T=Tp+ V2(AD) min 


Assim tem-se que na nova escala T', os conteúdos 
entálpicos das correntes quentes e frias podem ser 
livremente permutadas a um mesmo nível de 
temperatura; pelo que, em cada intervalo de 
temperatura, haverá um déficit ou um excesso de 
entalpia no processo, na sua totalidade, ditado pelo 
respectivo balanço entálpico nesse intervalo de 
temperatura. 


Como a entalpia flui ao longo do processo para 
temperaturas mais baixas, determina-se para cada 
temperatura extrema do intervalo o balanço 
acumulado de entalpia a essa temperatura por mera 
adição algébrica do balanço dos intervalos que lhe 
são superiores, como se apresenta na coluna mais à 
direita (ver fig.3). Este “o princípio que permite obter 
a Cascata de Calor do Processo!9). 


Este valor da temperatura (T') em que o processo tem 
maior deficiência da entalpia designa-se por "Ponto 
de Estragulamento "ou "Pinch" e, tem uma 
importância decisiva no traçado do diagrama de 
processos se, se quiser optimizá-lo sob o ponto de 
vista energético. Bastará dizer que toda a entalpia 
fornecida ao processo abaixo do ponto de 
estrangulamento é não só desperdiçada como 
contraproducente, bem como toda a entalpia retirada 
ao processo acima do “Pinch" ou "Ponto de 
Estrangulamento". 


A existência deste "Ponto de Estrangulamento” foi 
reconhecida teoricamente por Umeda et Alter, em 
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1978 (I0)que, dela não tiraram todas as 
consequências, tendo sido Bodo Linnhoff que, após a 
conclusão da sua tese de doutoramento no Dep. De 
Engenharia Química da Universidade de Leeds em 
1979, a lançou na realidade industrial da ICI como 
uma ferramenta indispensável da síntese do processo 
químico, com notáveis e brilhantes consequências 
práticas. 


5.2 Aplicação ao exemplo simples referido em 3. 


Para sabermos quais os consumos mínimos teóricos 
em utilidades quentes e frias, do processo referido no 
exemplo apresentado em 3., basta traçar a respectiva 
cascata de calor ; para isso basta assumir o valor 
indicativo do (AT) min dos permutadores de calor a 
instalar que por simplificação podemos continuar a 
assumir como sendo de 20ºC, 


CASCATA DE CALOR DO PROCESSO 


Tec) + Oq air DISXIO 


380 


a 
Dperin= MM0X 10h 


FIGURA 3 


Conclui-se portanto que o Consumo Mínimo Teórico 
da Utilidade quente deste processo químico é de 
335 X 103Kcal/h, a serem fornecidas acima de 250ºC, , 
dem vez do valor optimizado de forma heurística em 
3,1 de 660 x 103Kcal/h, €, 0 consumo Mínimo Teórico 
de utilidade fria é de 440x 103Kcal/h a serem 
fornecidas abaixo de 2500C, em vez das 
750x 103Kcal/h, obtidas em 3.3. 


Estes valores Teóricos são uma característica 
intrinseca do próprio Processo, baseados em 
princípios termodinâmicos e, não poderão ser 
alterados a não ser que o próprio processo seja 
alterado. Isto significa também que para que um 
Diagrama de Fluxos do processo esteja teoricamente 
optimizado não é possível fluir calor duma corrente 
acima do ponto de estrangulamento para uma 
corrente abaixo do ponto de estrangulamento. 


Para que o sistema esteja optimizado não é pois 
possível fluir calor da zona que absorve calor( acima 
do ponto de estrangulamento) para a zona que rejeita 
calor (abaixo do ponto de estrangulamento). 


Todo o calor X que o processo faça passar através do 
ponto de estrangulamento representa não só um 
excesso de calor que se tem que fornecer ao sistema 
acima do ponto de estrangulamento como um excesso 
de refrigeração que se tem de fornecer ao sistema 
abaixo do ponto de estrangulamento. 


Ou seja quanto mais calor for fornecido ao sistema 
acima do ponto de estrangulamento mais calor tem 
de ser retirado ao sistema abaixo do ponto de 
estrangulamento, o que representa uma dupla perda 
muitas vezes ainda agravada pelo aumento 
desnecessário do investimento nos permutadores de 
calor requeridos para electuar estas operações. 


Desde já se pode observar que no processo optimizado 
heuristicamente se está a transferir calor em P1, 
duma corrente acima do ponto de estrangulamento 
(390 - 2700C) para uma corrente abaixo do ponto de 
estrangulamento (50 - 1700C) o que claramente é 
uma situação não optimizada. 


5.3 Outras consequências práticas desta nova 
tecnologia. 


Para nos podermos guiar duma form mais precisa 
nas alterações concretas a efectuar no Diagrama do 
Processo com o fim de optimizarmos o consumo de 
utilidades vamos traçar o diagrama T - H do 
processo, considerando ainda como simplificação um 
(AT) min =20º€ , para os permutadores de calor em 
contracorrente a utilizar neste processo. 


Comecemos por recordar os diagramas elementares 
T-H, em que se representam correntes de calor 
quentes e frias, ou seja, aquelas que no processo têm 
de ser aquecidas ou arrefecidas. 


O CP representa a capacidade calorífica da corrente 
tal como já foi definido. Note-se que a inclinação da 
recta é inversamente proporcional ao CP 
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Se um mesmo processo Liver duas correntes quentes 
(ou frias) a curva que representa o processo integrado 
é simples de obter, dado que a variação da entalpia se 
pode definir em relação a qualquer estado de 
referência. 


Como é óbvio a construção do diagrama do processo 
integrado das correntes frias será idêntica. 

Para se construir o diagrama do processo integrado 
que inclua simultaneamente as correntes quentes e 
as correntes frias bastará representar as 
temperaturas na nova escala T", já definida atrás em 
5.1, que, permite fazer adições algébricas de 
conteúdos entálpicos de correntes a um mesmo nível 
de temperatura, conforme se verifica na Figura 4. 


Neste Diagrama vê-se claramente que todo o déficit 
da entalpia no intervalo 320-310º pode ser fornecido 
pelo excesso de entalpia gerado no intervalo 
380-3200C. Em termos processuais isso significa que 
a corrente [2] deverá fornecer calor à corrente [5] no 
permutador P.5, antes de ir aquecer a corrente [1] em 
Pa. 


REPRESENTAÇÃO GRÁFICA DAS NECESSIDADES ENTÁLPICAS DUM PROCESSO 


* Diagrama T = H duma corrente 


ar=qo 4 


abiagrama T- E de duas correntes quentes 


- Disgrame 1 = EH duma corrente quente em conjunto 


com uma corrente fria, numa escala T' 


FIGURA 4 


Obtem-se assim o Diagrama Processual Il com um 
consumo de utilidades quentes de 500x 103Kcal/h, a 
comparar com um Consumo Mínimo Teórico de 
335 x 103Kcal/h, 

Obtem-se assim a Curva Composta Global do 
Processo: 


Recordemo-nos que se estabeleceu como objectivo 
geral que não se deve instalar mais nenhum 
permutador de calor do que os existentes no 
Diagrama Processual I, ou seja, não deve ser 

aumentado o investimento. 


CURVA COMPOSTA CLOBAL DO PROCESSO 
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FIGURA 5 


- Poder-se-à, nestas condições, ir mais lunge na 
optimização do Consumo de Utilidades”?:- 


De [Poe 2s)s nino Lo! Eezalib 


- Relrigetagão total Fornecida pelas utilidades 


[Dbat, 47 L PO DO ja [LDA tás tosa LO! Ecalfi 


FIGURA 6 
Então atentemos agora no seguinte 


Se a corrente [2] após ter sido arrefecida em P.3 e Pl 
(e depois como é óbvio de sair de P.5) atingir 
exactamente os 180ºC isto significa que não necessita 
de ser refrigerada por utilidade fria e, que o 
permutador P.2 deixou de ser necessário. 
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Assim o permutador P.2 fica disponível para ser 
utilizado para fornecer o calor suplementar à 
corrente [3] a partir da utilidade quente de forma a 
que esta atinja os 260º€ . 


Temos assim o Diagrama Processual IV em que o 
consumo de utilidades quentes foi reduzido a 
410x 103Kcal/h e o consumo de utilidades frias a 
540 x 103Kcal/h e isto sem qualquer investimento 


suplementar. 


= E y ptnecidoa polida 
DIAGRAMA PROCITEVAL 17 ao E DSETIA PR — 


FIGURA 7 


Relativamente ao diagrama inicial alcançamos 
assim uma redução de 57,7 %, em termos de 
utilidades quentes. 


Mesmo em relação ao sistema optimizado 
heusticamente (Diagrama Processual 11) 
conseguimos alcançar uma poupança de 34,1 %, 
tmbém em termos de utilidades quentes 


O estudioso de engenharia química pode assim 
concluir que, com o mesmo equipamento básico ( um 
reactor e duas colunas ão) os resultados 

económicos do processo serão drásticamente 
influenciados pela forma como estes forem integrados 
ou sintetizados e, que a tecnologia do " ponto de 
estrangulamento” é de facto uma ferramenta 
indispensável na Síntese do Processo . 


Por outro lado e, uma vez que o custo dos 
permutadores é uma função da área de transferência 
de calor e que esta, por sua vez, é uma função da 
quantidade de calor a transferir, esta é mais reduzida 
no Diagrama Processual IV no seu conjunto, pelo que 
uma maior poupança de utilidades se [ez em 
silmultaneamente com um menor investimento, ao 
contrário do que até há bem pouco tempo era 
universalmente acei 


Consequências da tecnologia do Lo. 
estrangulamento" para o projecto de fornecimento de 
utilidades quentes e frias ao processo. 


Consumo de Utilidades 


Poupança 
j 
Quentes Frias Utilidade Utilidade 
Ucul/hi iRcalih) quente fria 


[Comeomeinimotesrco [385 [o 


(1) Relativamente so Diagrama Processuul | 


QUADRO 2 


A forma como se projecta o fornecimento de 
utilidades a um processo é por vezes decisivo para a 
respectiva viabilização económica. 


Isto é válido tanto em projectos feitos de raiz como, 
em projectos de modificação e optimização duma 
unidade já instalada. 


Se analisarmos de novo e atentamente a Curva 
Composta Global do Processo, em conjunto com a 
respectiva cascata de calor, verificamos que se 
optimizarmos a recuperação do calor entre as 
correntes do próprio processo, 0 vapor apenas precisa 
de ser fornecido a uma temperatura acima de 2800C; 
ao contrário, se não se tivesse feito essa recuperação ( 
ver Diagrama 1) o vapor teria de ser fornecido acima 
de 330º€, 

* E, se o vapor disponível tiver uma tem 
superior a 280007 


tratura 


Então , caso tal seja possível. esse vapor deverá 
accionar primeiro uma turbina de contrapressão, de 
forma a que o respectivo vapor exausto esteja em 
condições de fornecer calor ao ebulidor da Coluna de 
Destilação 1. 


*E qual a utilidade fria a utilizar ? 

Como se vê no Diagrama IV a utilidade fria apenas é 
utilizada nos condensadores das Colunas de 
Destilação 1 e 2, com temperaturas finais de, 


respectivamente, 120 e 1800€, 


A água de refrigeração é pois uma solução correcta e 
provavelmemte óbvia. 


5.4 Perspectivas futuras 


+ 


Integração de Processos quanto à 
duma bomba de calor num processo ? 


Convirá antes de mais clarificar à luz da tecnologia 
do “ponto de estrangulamento”, qua! deverã ser a 
localização duma bomba de calor. 
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so — 


Como sabemos uma bomba de calor retira entalpia ao 
processo a uma temperatura Ta e devolve-a a uma 
temperatura Tg>TaA, consumindo para o efeito 
energia mecânica. 


Dado que através de meios exteriores ao processo se 
deve apenas fornecer calor acima do ponto de 
estrangulamento e apenas retirar calor abaixo do 
ponto de estrangulamento, então Tp>Tpe. 
Ta<Tpe. 


Uma localização apropriada para uma bomba de 
calor poderia ser, aparentemente, entre a corrente [6] 
a 1800C (T') e a corrente [5] a 3109€ (T'”). 


* Todavia, será esta localização da bomba de calor 
uma solução correcta? 


Analisemos de novo a Curva Composta Global do 
Processo : 


- À corrente [5] pode ser aquecida com a corrente [2] 
como já se verificou na prática no Diagrama 
Processual IV; 


- e a entalpia disponível na corrente [6] pode ser 
utilizada para pré-aquecer a corrente [1] (como se 
verá caso se pretendesse atingir o Consumo Minimo 
Teórico de Utilidades). 


- Conclui-se assim que, embora à partida não seja 
errada, a colocação desta bomba de calor não faz 
sentido em termos de optimização do consumo de 
utilidades . 

De facto a análise da Curva Composta Global permite 
estabelecer duas condições que têm que ser 
satisfeitas em simultâneo para justificar em termos 
económicos a localização duma bomba de calor: 


a corrente fria, situada acima do ponto de 
estrangulamento, terá que receber entalpia duma 
utilidade; 


- à corrente quente, situada abaixo do ponto de 
estrangulamento, terá que ser refrigerada 
obrigatóriamente por uma utilidade. 


* Então onde instalar uma bomba de calor, no 
processo que temos vindo a analizar? 


Por exemplo retirando calor à corrente |4] para 
fornecer à corrente [3] Tratar-se-à pois duma bomba 
de calor que retira ao condensador para o fornecer ao 
ebulidor da Coluna da Destilação 1 e que, conforme se 
pode ver no Diagrama da Curva Composta Global, 
substitui utilidades. 


ue outros aproveitamentos potenciais de calor 


Se analisarmos cuidadosamente a Curva Composta 
Global verificamos que a corrente [2], no Diagrama 
Processual IV, deve abandonar P.3 a 2600ºC e, não a 
3000ºC, para se conseguir fazer a Recuperação 
Máxima Teórica. 


Dessa forma a entalpia remanescente na corrente |2] 
é insuficiente para pré-aquecer a corrente [1] Será 
assim necessário utilizar o vapor? 


Verificamos fácilmente que não, será necessário sim 
utilizar para isso parte da entalpia da corrente [6]. 


6. CONCLUSÃO 


O desenvolvimento de técnicas de Sintese eficientes 
são uma parte essencial do curriculum de qualquer 
curso de engenharia quimica. 


Exemplos como o desenvolvido neste trabalho são 
simultaneamente simples e estimulantes, visto que 
eles usam algumas das tecnologias mais recentes e, 
constituem uma forma muito importante para o 
estudante usar a sua imaginação duma forma 
criativa, com respeito absoluto pelos conhecimentos 
que ele terá de ter relativamente às uvperações 
unitárias envolvidas no processo. 


Dai a sua inclusão no programa da cadeira de 
Optimização de Processos, leccionada a partir da 
recente alteração curricular aos alunos do 2º 
semestre do 4ºano, do curso de engenharia Química, 
ramo de Processos e Indústria, do Instituto Superior 
Técnico. 
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ABSTRACT 


It is shown that when the upper-shielding of a 
conductor-backed dielectric slab is backed to infinity, 
only the surface modes of the corresponding open 
condocter-backed slab with the same dielectric 
thickness are obtained, if and only if the original 
closed waveguide modes are slow-modes. The 
continuous spectrum is not obtained in this limit, 
therefore making the so-called discretizing 
procedure an approach to be used with precautions. 


INTRODUCTION 


Open dielectric waveguides are generally analyzed 
by the usual procedure of discretizing the continuous 
spectrum by closing the original waveguide with 
perfectly-conducting walls [1-3]. The advantage is 
obvious: for closed waveguides complete field 
representations require only discrete modes [4]. 
This is specially convenient when a mode-matching 
technique is used for the analysis of transversely 
discontinuous dielectric waveguides [5]. 

In this short paper, the modal equations of the 
condoctor-backed  dielectric slab with 
upper-shielding (Fig. 1) and of a similar structure 
but without the upper-shielding (Fig. 2) are 
compared .“The bottom line in the present work is to 
analyze the case of the upper-shielded structure of 
Fig. 1, when the perfectly-conducting wall placed at 
x=a+t is backed to infinity, i.e. when a >. With 
this simple example, the consequences of the usual 
discretizing procedure can be analyzed... 


COS Rss es Es 


Original recebido para publicação em 17.02.88. 


RESUMO 


Mostra-se que quando se faz recuar até ao infinito a 
parede superior de um guia de planos paralelos 
parcialmente preenchido, apenas se obtêm os modos 
superficiais da correspondente placa dieléctrica 
assente em plano de terra e com a mesma espessura, 
se e só se os modos de guia fechado original são 
lentos. Neste limite, o espectro contínuo não é obtido, 
tornando assim o chamado método de discretização 
uma técnica a ser utilizada com cuidado. 


Fig. 1 - Conductor-Backed Dielectric Slab With 
Upper-Shielding. 


DISCUSSION 


For the open waveguide of Fig. 2 let h and a be, 
respectively, the transverse wavenumber in medium 
(1) and the transverse attenuation in the air 
(medium (2) ). Hence 


h — ko(Er- Eefl qua 


a= k., Cn =)" (1) 
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where £r is the relative dielectric constant of medium 
(D), ko= wleouo)” and cem is a relative effective 
dielectric permittivity of the guided mode which is 
related to the longitudinal wavenumber k by 


tetr= k2/k,2 (2) 
x 
PR o (2) 
k A (1) 


—e ge 


Fig. 2 - Conductor-Backed Dielectric Slab. 
Defining 


u= hl 
w= at (3) 
v= koble,-1)12 


where t is the thickness of medium (1), one has 

tu + w2 =v2 (4) 
both for TE and TM surface modes. 
Moreover, one has 

w= - ucotgu (5) 


for TE surface modes and 


| 
w=— utgu 


E 
x 


(6) 
for TM surface modes [4]. One should note that w and 
u are always positive or null, 1.e. |I=eam <te,. Due to 
the boundary condition at x=0 (E,=E,= 0) only 
surface modes with H, (E,) even (odd) functions of x 
for TM (TE) modes can propagate in this open 
waveguide. Therefore, there are only TM, and TE, 
surface modes with p even and q odd. For all surface 
modes the cutoff wavelength A corresponds to w= 0. 
Hence, for TM, surface modes, one has 


(7) 


with p even. For p=0, Eq. (7) leads to (Aç) 9=, In 
fact, the TM, is the only surface mode without cutoff. 
For TE, surface modes, with q odd, (Ac) Is still given 
by Eq. (7) ifone replaces p for q. 


However, a complete field representation is not 
possible in terms of these surface modes [4]. 

For the closed waveguide of Fig. 1, Eq. (4) is still 
observed |5]. 

One may write [5], for TE modes, instead of Eq. (5), 


a 
ucotg vu + w coth( r, w)= 0 


(8) 
and for TM modes, instead of Eg. (6), 
l a 
-ute u-w tanh( à w)= O 
£ 
i (9) 


One should note that for the closed-waveguide of 
Fig.1, w may be positive (£egr > 1) or imaginary (ger 
<1). In faet, for TM, modes (with p even) one obtains 
from Eq. (9) that geg> 1 for 


t p 


> ps 
À at 1)? 


(10) 


As for TE, modes (with q odd), according to Eg. (8), 
tem >1 for VA> (UM min where 


ante 1)? E À 

(LT) 

Hence,according to Eq. (10), the TM, mode is a 
slow-mode,i.e. its phase velocity 


(12) 
is lawer than c= 1/(e,H9)!2 for all values of t/A. Apart 
from the TM, mode,all other modes are fast (e.g <1) 
for low values of t/A and slow ( cer >1) for high 
values of tA. | 
One should also note that, as for the open waveguide, 
the TM, mode of the closed waveguide has no cutoff. 
In fact, according to Eq. (9), for a/t constant, one has 


e la/t+ 1) 
li 


a A 


———— & 
- e att] 

(13) 
where the equality stands for e, = 1. When a > this 
limit becomes 


lime,y = 1 
a —>00 (14) 
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which is also the limit for geg when A —0 for the 
TM, mode of the open waveguide of Fig. 2. Thus, the 
TM, mode is the fundamental mode for both 
structures of Figs. 1 and 2. For higher order modes, 
there is always a cutoff wavelength Aças for the open 
waveguide. However, it corresponds to eeg=0 (k=0) 
and is given by 


(15) 
both for TE and TM modes, according to Egs. (8) and 
(9). Hence,for eeg=0, TMp and TE,.1, with peven and 
not null, are degenerate. 

If one considers a slow-mode of the closed waveguide 
of Fig. 1, according to Eqgs. (10) or (11), one has w>0. 
Hence,one can see that when a >, the limits of Egs. 
(8) and (9) are, respectively, Eqgs. (5) and (6). 
Therefore, if t/A observes Egs. (10) or (11). ie. if a 
given mode of the closed waveguide is a slow-mode, it 
is reduce to Lhe corresponding surface mode of the 
open waveguide when a —0. 

However, for the closed waveguide of Fig. 1, the 
considered TE and TM modes form a complete field 
representation as the field-wave equation together 
with the boundary conditions constitute a 
Sturm-Liouville system [4]. Therefore, when a —0 
and under appropriate conditions, what was a 
complete field representation for the close waveguide 
originates an incomplete field representation for the 
open waveguide. 


CONCLUSIONS 


It has been shown, by means of a simple example, 
that the so-called discretizing procedure of the 
continuous spectrum of open dielectric waveguides is 
an incomplete method for the analysis of these 
waveguides. The example studied in this short paper 
clearly reveals that only the discrete eigenvalue 
spectrum of the open waveguide of Fig. 2 may be 
obtained by this procedure and only under 
appropriate conditions. The continuous eigenvalue 
spectrum is lost with the discretizing procedure. 
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ABSTRACT 


A computer based technique to measure the 
correlation of velocity signals from a laser 
anemometer and temperature signals from a 
compensated, small-diameter thermocouple is 
described together with the use of the technique in 
premixed flames. It is shown that the probability 
density function of measurements conditioned by the 
arrival of velocity signals is close to 
density-weighted. 


1. INTRODUCTION 


The purpose of the work reported here is to describe a 
method for measuring the correlation between 
velocity and temperature fluctuations in a turbulent, 
premixed flame using laser-Doppler velocimetry and 
digitally -compensated bare-wire thermocouples. 
There are several incentives for making such 
measurements. For example, the correlation 
between velocity and temperature fluctuations is 
closely related to the turbulent transport of the 


E a dana and A 
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products of combustion, which has been shown to 
occur in the "countergradient" sense by recent 
analytical work of Bray et al (1981), Masuya and 
Libby (1981) and Libby et al (1985), but detailed 
comparisons with experiments are required. Also, 
frequency resolved values of velocity and 
temperature or density will allow conclusions to be 
drawn as to whether Favre-, or Reynolds-averaged 
data is measured (e.g. Bilger, .1981). 

The instantaneous measurement of velocity can be 
made by laser-Doppler velocimetry. It involves the 
measurement of the frequency of light scattered from 
small particles added to the flow: the output of the 
laser velocimeter is, therefore, not continous and 
depends on the presence of a particle in the control 
volume formed by the laser beams. A reliable 
measurement of temperature cannot be said to exist 
and the method chosen here is to use a fine-wire 
thermocouple. 

The thermal inertia of the thermocouple gives rise to 
first-order damping of the frequency response, for 
which compensation must be made. The time 
response of a bare-wire themocouple can be described 
by the following diferrential equation. 
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dT, 

C— + T, = aa 
Where Tp is the instantaneous temperature of the 
gas, Th is the intantaneous measured bead 
temperature and / the "time constant” of the 
thermocouple. There is an axtensive body of 
literature, for example Ballantyne an Moss (1977), 
Yule, Taylor and Chigier (1978),Yoshiba and Tsuji 
(1979), Lenz and Guenther (1980), Lockwood and 
Moneib (1980), Yanagi and Mimura (1981), Tanaka 
and Yanagi (1983), Sivasegaram and Whitelaw 
(1983),LaRue, Samuelsen and Seiler (1984) and Nina 
and Pita (1985), describing how such compensation 
can be affected, once the time constant of the 
thermocouple has been found. Since this constant 
depends on the heat transfer coelficient, its value is a 
function of the local velocity, temperature and 
composition. 
Our objective is to describe a computer data 
acquisition system to obtain joint probability density 
functions of velocity and temperature, their 
correlations and conditationally-sampled values. 
Sample results are presented to help investigate the 
kind of averaging process perfomed by the 
anemometer. A detailed description of the technique, 
which includes the effect of velocity and temperature 
on the time constant of thermocouple, has been 
provided by Heitor (1985) and Heitor et al (1985). 
The following section describes the experimental 
method and the arrangements for the digital 
recording of data and its subsequent reduction by 
laboratory computers. The third section presents the 
results of investigations to evaluate the technique. 
The paper concludes with a summary of the main 
findings. 


2. EXPERIMENTAL METHOD 
2.1 FLOW CONFIGURATION AND 
INSTRUMENTATION 


The flame used in the present study was stabilised on 
a disc of 0.056m diameter, D,which was positioned at 
the exhaust of a 1.83m long straight pipe of 0.080m 
diameter. The reactants were air and natural gas 
and the annular bulk velocity at the trailing edge of 
the disc was 18 m/s, resulting in Reynolds numbers 
based on the pipe diameter of about 60 x 103: the 
equivalence ratio of the mixture was maintained at 
about 0.8. The air flow was seeded with aluminium 
oxide particles from a spark discharge particle 
generator, in order to provide light scattering centres 
for a laser-Doppler velocimeter. 

Velocity was measured by a dual-beam, 
foward-scatter arrangement with sensitivity to the 
flow direction provided by light-frequency shifting 
from acousto-optic modulation at 40 MHz. The 


scattered light was collected by a lens and focussed 
onto the pinhole aperture of a photomultiplier, the 
output of which was mixed to give a net frequency up 
to 9 MHz. The signal was amplified (20 dB) and 
measured by a frequency counter. The output of the 
counter was passed (via an IEEE - 488 interface bus) 
to a 16-bit laboratory mini-computer with 32K 
storage words. 

Temperature was measure by thermocouple 
junctions made by butt-welding platinum wires with 
platinum-rhodium (13%) wires of either 0.015 or 
0.040mm diameter. The output of the thermocouple 
was differentially amplified (x100) by a low noise, 
low drift amplifier and digitised by a 12 bit 
analogue-to-digital convert at sampling rates of up to 
about 38 KHz. The noise level could be kept below 
ômv (0.098% of full scale deflection), corresponding 
to a maximum temperature error of +t2 K. The 
samples were stored in a memory of an eight-bit 
laboratory microcomputer whith 48 k storage words. 
The time constant of the thermocouple was 
evaluated numerically and on-line, as well as the 
subsquent compensation for thermal inertia. This 
have included the dependence of the time constant on 
the velocity and temperature of the flow. 


2.2 SIMULTANEOUS DATA LOGGING 


The algotithm for compensation requires (see below) 
a temperature record of at least 0.5ms duration, 
sampled at above 17 KHz before and afther the 
occurrence of a velocity measurement..This 
requirement together with available computer 
resources, precludes the continuous recording of 
temperature over the time for occurence of a 
statistically acceptable number of velocity 
measurements, say over 1000. The most economical 
use of the limited memory space of the computer 
which is dedicated to the measurement of 
temperature is to store the samples only when a valid 
velocity measurement is made. The following 
paragraphs describe how this was .accomplished 
together with the flow chart of figure 1. 

The two computers were linked by two separate lines 
which were used to synchronize the measurement of 
velocity and temperature and to pass the 
temperature data to the 16-bit machine for 
subsequent reduction and presentation. One line was 
a standart (EIA-RS232C) serial line, which allowed 
communication in either direction, and the other was 
an analogue line which was controlled by the larger 
machine and read by the eight-bit machine. 
Separate, and initially independent, programs were 
executed in each-machine which initialised the 
various interfaces and prepared for the recording of 
data. The serial communications interface was used 
to ensure that both computers had completed this 
phase before any measurements were taken. 
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The recording of data was contingent on the presence 
of a valid velocity measurement. In the absence of a 
measurement from the frequency counter, the 16-bit 
computer was idle while the smaller machine 
continuously sampled the output of the thermocouple 
and stored the result in a circular buffer memory (the 
voltage output of the thermocouple was converted to 
temperature by means of a "look up” table in the 
8-bit computer).The occurence of a velocity 
measurement was signalled by the counter to the 
large computer (via the IEEE-488 interface bus) and 
the latter immediately performed the following two 
operations. First, a pulse was issued on the analogue 
line (see above) as a signal to the 8-bit machine that 
a velocity measurement had been acquired. Next, the 
16-bit machine checked the number of zero crossings 
and the gate time, measured by the counter, to 
determine whether these satisfied the validation 
requirements. The outcome of the check was 
subsequently signalled to the 8-bit machine (via the 
RS-232 line). A valid velocity measurement resulted 
in the transfer of the 20 most recent temperature 
samples from the circular buffer to a separate 
storage area of the 8-bit computer and in the storage 
of the Doppler frequency in the 16-bit computer. This 
procedure bracketed each velocity measurement 


with ten temperature samples before an after the 
occurence of velocity. When both these memory 
operations had been completed, the cycle was 
repeated until 512 valid data points had been 
obtained, and the temperature samples were then 
transferred from the smaller to the larger computer. 
A compiled FORTRAN program compensated the 
temperature measurements with the compensated 
temperature, T., given by a 


dt, 


P=T+ TU) 7 


with the temporal derivative calculated by central 
differencing from the temperature samples which 
bracket the velocity measurement (see above). It is 
an imposed restriction that the derivative be 
calculated between samples which have a 
temperature difference, at least, between 1 and 4 
kelvins so that the derivative does not reflect the 
effects of single bit changes of the least significant 
bit of analogue-to-digital converter. This imposition 
requires the recording of more than just two 
temperature samples (one just before, and one just 
after, the occurence of a velocity measurement) and 
in pratice twenty samples were collected (see above). 
The process was repeated until a statistically 
significant population had been achieved and the 
moments of velocity and temperature were then 
calculated and plotted. 

This computer data acquisition system limited the 
conditional sampling frequency to about 73 
measurements/sec, which is larger than the 
maximum data rate recorded in the flame (about 50 
measurement/sec) bedfore a significant deterioration 
of the thermocouple bead by seeding-particle 
accretion occurred. 


3. SOME RESULTS AND DISCUSSION 
ONTHE AVERAGING PROCESS. 


The foregoing described the arrangement of the 
instrumentation which was used to measure velocity 
and temperature simultaneously. The purpose of this 
section is to investigate the king of averaging process 
which is performed by the anemometer. Heitor et al 
(1985) have compared probability density functions 
of continuous records of temperature with those 
conditioned by the occurence of valid velocity 
measurements and showed that the latter may be 
considered close to Favre-weighted without the need 
for further densi ty weighting (1). The main 
approximation involved in this statement is that no 
account has been taken of the relation between the 
probability of the occurrence of the velocity signal 
and the magnitude of velocity ("velocity biassing", 
e.g. Dimotakis, 1976, Edwards and Meyers, 1984). In 
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other words, the Favre-weighted probability density 
functions are biased but for the turbulence 
intensities considered here, the errors are less than 
+4 and -4% for the mean and variance respectively 
(Glass and Bilger, 1978). 


Figure 2 shows axial profiles of the mean and 
normalised variance of progress variable along a line 
within the shear layer of the disc stabilised flame 
described previously. The values were obained from 
four averaging processes as follows: (1) a continuous 
record of temperature, which gives rise to an 
unweighted average of temperature; (ii) 


density-weighted values of the continuous record of 
temperature; (iii) temperature measurements 
conditioned by the arrival of valid velocity signals at 
the same point; and (iv) the density-weighted values 
of the conditioned temperature measurements. The 
results support the conclusions of Heitor et al (1985) 
based on the analysis of probability density 
distributions: measurements conditioned by the 
existence of a valid velocity sample correspond to 
“automatic” Favre averaging due to the different 
densities of the seed particles in the hot products 
comparated to the relatively cold reactants. The 
results allow to quantify the differences between 
Favre, or density-weighted, and unweighted 
temperature moments: the maximum diferences for 
the first moment lie in the reaction zone with values 
up to about 10% of the local Favre-averaged 
temperature, for intensities of the temperature 
fluctuations around 25%. The associated differences 
in the r.m.s. of temperature fluctuations are about 
20% of the local-averaged r.m.s. valu 


4. SUMMARY 


This work evaluated real-time algorithms which 
have a general form to be applicable in different 
measurement environments. The explicit forms of 
these algorithms are presented and discussed with 


emphasis on the measurement of unconditioned and 
conditioned velocities and correlations of velocity 
and temperature using laser-Doppler velocimetry 
and bare-wire thermocouples in premixed flames. An 
important finding of the work reported in this paper 
is that the probability density function of 
measurements conditioned by velocity is close to 
density-weighted. 

In general, the method presented facilitate on line 
conditional sampling in applications where 
conditioned variables are used to carry out and 
validate the measurements of other quantities. 
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RESUMO 


O poliacetileno é um polímero orgânico cuja 
condutividade eléctrica pode ser alterada, quer 
química quer electroquimicamente, com dopantes 
que aceitam e doam electrões, conferindo-lhe 
propriedades características dos semicondutores e 
dos metais. Esta propriedade despertou um enorme 
interesse na comunidade científico-tecnológica face 
às potencialidades práticas decorrentes, como se 
descreve nesta comunicação. 


1. INTRODUCTION 


A new chapter in the evolution of conducting 
polymers began with the 1973/1975 discoveries of 
metallic and superconducting properties of polymeric 
(SN); [1]. When in 1977, MacDiarmid, Heeger, 
Shirakawa,and coworkers [2,3] demonstrated the 
semiconducting and metallic properties of the linear 
conjugated organic polymer polyacetylene, a surge of 
activity began to synthesize and characterize 
members of this new class of organic electronic 
materials. 
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ABSTRACT 


Polyacetylene is a typical semiconducting and 
metallic organic polymer,whose electrical 
conductivity can be tailored with a variety of electron 
accepting or electron donating reagents. This 
unusual characteristic has opened an extensive area 
of research on its possible technological applications, 
as described in this paper. 


Since then other conducting polymers such as poly 
(paraphenylene)[4], polypyrrole [5], polythiophene 
[6], and polyaniline [7] have been discovered which 
have also attracted the increasing attention of the 
whole scientific community, because of their unusual 
electrical conductivity and possible technological 
applications, which include reversible electrodes for 
batteries, electromagnetic interference shielding, gas 
sensors, displays, conducting plastics, capacitors, 
junction devices and metals replacement. 

This article describes briefly the preparation and 
characterization of thin films of polyacetylene and 
reviews many ofits potential applications. Particular 
emphasis is placed on its electrochemical 
applications. 
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It is noted that by far the greatest amount of 
published material on conducting polymers has dealt 
with polyacetylene films. 


2 PREPARATION AND CHARACTERI-. 
ZATION 


Polyacetylene is the simplest linear synthetic 
polymer. Fundamental unit consists from one carbon 
and one hydrogen atoms, as has been called simply 
(CH), Each carbon is o bonded by sp? hybrid orbital 
to one hydrogen and two adjcent carbon atoms to 
form a planar chain molecule. The remaining fourth 
valence electron of the carbon atom is in P, orbital, n 
electron, which lies perpendicular to the molecular 
plane. Thus, the molecule forms a 
quasi-one-dimensional lattice of the nm electrons. 
Such polymer has been the subject of many 
theoretical investigators [8-15] on the electronic 
structure and bond alternation. Depending on the 
degree of the bond alternation,polyacetylene was 
expected to be an insulator to semiconductor in its 
electrical conductivity [16]. 


The technique for polymerisation developed by the 
Japanese employs a so-called Ziegler catalyst - a 
mixture of Ti (0C4H9)4 and Al (C9Hs5)s, first reported 
by Giulio Natta and co-workers at Milan Polytechnic 
in 1958. The Japanese groups directed acetylene on 
to the surface of the catalyst solution where a black, 
shiny film of polyacetylene from 10-5 to 0.5 em thick 
formed. An alternative method, which can be used to 
make films with a large area, is to spread a thin 
layer of catalyst solution over the walls of the 
reaction vessel. When acetylene gas enters the vessel 
a polyacetylene film forms on the wet surface. The 
polyacetylene produced in this way has a smooth, 
shiny surface in contact with the reaction vessel and 
a rough, dull surface on the interior. Although 
usually referred to as a film,the polyacetylene is a 
low-density sponge-like structure made out of fibrils 
with diameters of a few tens of nanometres. 


Polyacetylene has two isomeric forms, trans and cis, 
or as a mixture of the two forms depending on the 
preparative conditions. Thus, if polymerization is 
carried out at Dry Ice temperature (ca. -78ºC), the 
film is formed almost completely as the cis isomer, on 
the other hand, ifa temperature higher than 100ºC is 
used, the film is formed as the trans isomer. As 
shown in Fig. 1, the room temperature DC 
conductivity and the activation energy for 
conduction of the films depend upon the cis-trans 
contents, varying from 10-9 S/em and 0.5 eV for the 
cis isomer, and 105 Sem and 0.3 eV Tor the trans 
isomer respectively. 

Polyacetylene, either cis or trans isomer, can be 
doped chemically or electrochemically with a variety 


of electron accepting reagents and electron donating 
reagents to form p- and n-type semiconductores, 
[CH”'(Adyly and [CH*(D'),]y where A end D 
represent acceptor and donor dopant species, 
respectively, and y is the dopant concentration [1, 
17]. The conductivity of doped polyacetylene can be 
enhanced dramatically, and it is this property that 
stimulated interest in this material. For example, 
additional electron-accepting molecules, such as the 
halogens bromine, chlorine, iodine, or arsenic 
pentafluoride, caused the conductivity of 
trans-polyacetylene to increase from 10-5 S/cm to 
over 500 S/cm, the final value being close to that of 
metals, as Fig. 1 shows. Conductivities higher than 
those of some metals like iron and platinum. have 
recently been reported by Naarmann's group, of 
BASF [18]. The transport mechanism in the doped 
polyacetylene has remained controversial [1, 18-21]. 
It has been proposed that there are at least three 
dopant regimes [22]: 


“—. Semiconductora Se Metals — 


—— Insulators —s 


Drganie 


compounds é E (CH), ás á (CH) z 


Inorganic 
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Elemento ! ailasnd 51 Ge g Fe Cu 


Conduetivity (5 sa) 


Fig. 1 - The conductivity of (CH), compared with 
values for other materials. 


(1) y € 0.0001, the transport is via carrier hopping, 
mobility is small and activated; (ii) 0.001 ( y (0.01, 
the transitional region; and (iii) y ) 0.01, the metallic 
state. 


3. POTENTIAL APPLICATIONS 


The review presented in this section should serve to 
indicate the usefulness of doped polyacetylene films 
in the fields of chemistry and electrochemistry, 
energy, electronics and microelectronics,aeronautics, 
pharmacology and biology, and molecular 
engineering. 


3.1 Semiconductor devices 


Heterojunctions fabricated by polymerization of a 
thin (1-10 u m) (CH), film directly onto a 
single-crystal n-CdS substrate, and by pressure 
contact of a free-standing (CH), film (= 100 yu m) 
onto a single-crystal n-CdS surface, have been 
fabricated and characterized by current versus 
voltage (I-V), capacitance versus voltage (C-V ), and 
photovoltaic action spectra [23]. Rectifying 
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behaviour characteristic of a p-n junction diode has 
been shown by the I-vs-V data, thus confirming that 
as-grown polyacetylene is p type. From C-V 
characteristics, a residual acceptor concentration in 
as-grown films of p(CH), of 2x1018 cm3 has been 
inferred, and it has been concluded that in spite of 
the complex fibril morphology, (CH); can be treated 
as an effective homogeneous medium. The existence 
of a well-defined deep trapping state in polyacetylene 
with an energy near the center of the gap, has been 
implied by detailed studies of the photovoltaic 
response at energies below the energy gap for (CH),. 

Other heterojunctions and Schottky barrier diodes of 
polyacetilene with either n-type inorganic 
semiconductors or low work function metals have 
been reported. Homojunctions of p-type (CH); /n-type 
(CH), have also shown good rectifying 
characteristics [24-27]x. 

Application of polyacetylene to junction devices may 
be limited by the morphology of polyacetilene. 
Because of its fibrillar morphology it is very difficult 
to prepare very thin films without pin holes. In 
addition to this difficulty, fibrillar morphology 
prevents an effective contact with substrate. 


3.2 Solar Cells 


While present materials with potential application 
to solar cells have far less efficiency than the 10% 
now possible with silicon solar cells, the scientists 
are encouraged that it may be possible to produce 
lightweight, low-cost solar cells of polyacetylene. 
This is generated by the fact that it is easy to prepare 
thin polyacetylene films of large area with low cost, 
the absorption spectrum of the polymer is close to the 
solar spectrum, its absorption coefficient is very 
large (= 109 em-!), and its optical band gap is about 
l.5eVI1). 

A p(CH),: n-Z2nS heterojunction solar cell has been 
fabricated by Ozaki et al. [28] with open-circuit 
photovoltage of 0.8V. The utilization of(CH),, as the 
active photoelectrode in a photoelectrochemical cell 
using a sodium polysulphide solution as electrolyte, 
has been reported by Chen et al. [29].. A significant 
photoresponse has been demonstrated with 
open-circuit voltage Voc = 0.3 eV, and short-circuit 
current Is: =40 uA/cm?2, A schematic diagram of the 
PEC cell configuration is shown in Fig. 2(a) under 
ilumination of approximately 1 sun. Isopropyl 
alcohol has been added to the aqueous solution to 
achieve improved contact at the hidrophobic (CH); 
surface. À schematic diagram of the PEC cell 
configuration is shown in Fig. 2(a). Fig.2(b) shows 
the electrode processes occurring in the cell. Since 
(CH), is p-type, photogenerated electron-hole pairs 
are separated at the (CH), -electrolyte interface with 
electrons injected into the electrolyte. The resulting 
S2- irons migrate to the Pt eleetrode where they loose 


electrons to reform the S,2--species. The liberated 
electrons return to the (CH), electrode via the 
external circuit. It was shown that the polysulphide 
solution acts as an effective redox couple, being 
anticipated that a variety of redox systems should 
work in combination with the (CH), electrode. 


p- (CH) x PE 


modes 27 e 9 4-2 


Cathode: É + le” — 25º 
(a) (b) 


Fig.2 - Schematic diagram of PEC cell configuration 
(a) and electrode processes occurring in the 
cell (b). 


Other solar devices have been described in the 
literature [30, 31] but their characteristics are not as 
good as those of silicon p/n junction solar cells. 


3.3 -Batteries 


The fibrillar morphology & large surface area of 
(CH), prepared by the Shirakawa technique is 
presently believed to play an important role in its 
potential use as an electrode-active material. Since 
electrochemical reactions occur at the 
polymer-electrolyte interface, the larger this 
interface per unit mass, the greater should be the 
rate at which a (CH), electrode can be charged and 
discharged. The large void space between fibrils in 
(CH), films allows for ready access of electrolyte to 
(CH) in the interior ofa film, 

Free-standing 0.1 mm cis-(CH), films have been 
most commonly used in electrochemical and 
battery-type studies. The film is usually placed on a 
metal mesh gauze of platinum (current collector) 
which is then folded tightly upon itself so as to 
encase the film. The assembly then serves as the 
electrode. Cis-(CH), is usually preferred in 
constructing an electrode, since it is more flexible 
than the (more thermodynamically stable) 
trans-isomer. Since the cis-isomer is converted to the 
trans-isomer on electrochemical p- or n-doping[32)., 
the cis-isomer used in fabricating an electrode will be 
converted to the trans-isomer after one or two 
electrochemical cycling studies. 

When the positive and negative terminals of a DC 
power source are applied to a (CH), and a Li 
electrode, respectively, in an electrolyte such as a 
solution of LiCl04 in propylene carbonate (PC), 
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electrochemical oxidation proceeds readily. The 
electrode reactions are 
(CH); — (CHy +), + xyer 


at the (CH), (anode) and 
xyLit + xye — xyLi 


at the Li (cathode), resulting in the generalized 
overall charging reaction 


(CH), + xyLit — (CH”"), + xyLi 


For a specific anion, for example, Cl0-4y |, the 
resulting net reaction is 


(CH); + xyLi+CIO-4 > [CH”*) (CIO qy1y + xyLi 


Ifthe power supply is removed and the two electrodes 
are connected by a wire, spontaneous 
electrochemica!l reactions (discharge 
reactions),which are the reverse of those given above 
occur as electrons pass through the wire from the Li 
to the [(CH**) (Cl0-4)],, regenerating (CH), and 
LiCl04. These types of reactions provide the basis for 
the use of [(CHº"JA, | as a possible cathode-active 
material in rechargeable batteries. Cells of this type 
employing electrolytes consisting of salts of Cl0 4, 
AsF"6, BF'4, CI and CF35S0'3 dissolved in propylene 
carbonate or sulfolane have been studied [33-38 ] 
Solid-state batteries employing oxidized 
polyacetylene as the cathode-active material have 
also been reported [39, 40]. The characteristics of 
these cells are not particularly ideal, but they may be 
useful for powering electronic devices and liquid 
crystal displays. 


(CH), can also be reduced directly in cells of the type 
Li/Electrolyte/(CH);; then n-doped (CH), containing 
Lit dopant cations can be obtained: 


(CH), + xyer — (CHy'); (cathode) 


xyLi — xyLit + xye- (anode) 


(CH); + xyLi — [Lit (CH”)], (overall reaction) 


The n-doped (CH), can act as an anode in a battery 
when, durind the discharge reaction of a cell, it 
becomes oxidized, that is, the negative charge 
residing on the CH unit is reduced in magnitude. À 
simple example of a cell in which polyacetylene acts 
as an anode is that using a [Lit (CHy-)|, anode and a 
TiSo cathode in à LiCIQ (THF) electrolyte [41]. The 


net discharge reaction is therefore: 


[Li ACHY) + (TiS9)y > [Lit y(TiSoy)])z + (CH), 


and this cell has shown many desirable features of 
the doped anode as compared with the metal anodes 
(ie. Li). 

IH appropriate energy densities and recyclability 
could be attained, a most attractive rechargeable 
battery would be one in which both anode and 
cathode were fabricated from (CH), [42]., but a far 
more interesting cell (ie. with higher energy density) 
involves the use of a (CH”), anode and a (CH”"), 
cathode. In a typical example, in a cell consisting of 
two polyacetylene (CH), electrodes immersed in a 
solution of (n-Bu4N +) (PF 6) and allowed to charge to 
7% doping, the following electrochemical reactions 
are involved |43): 


[(n-BusN +)o o7(CH- 00], — (CH); + 


0.07x(n-BusN +) + 0.07xe- (anode) 


[(CH+0.07) (PF go 07)x + 0.07xe — (CH); + 


0.07x(PF-6) (cathode) 


[(n-BuaN +)o 07(CH-0.07)], + [(CH+0.07) 
(PF go 07x — 2 (CH); + 0.07(n-BugN +)(PF-6) 
(overall reaction) 


All-plastic batteries with polyacetylene electrodes 
have also been described [44, 45]. Plastic batteries 
using composite electrodes(obtained for example by 
polymerization of thiophenes on undoped 
polyacetylene) seem to possess higher efficiencies 
than those using polyacetylene or doped 
polyacetylene electrodes [46]. 


The study of batteries employing polyacetylene 
electrodes is still in its infancy. Consequently 
polyacetylene batteries are not already in 
commercial production, but it seems likely that 
commercial versions of batteries with polyacetylene 
electrodes will appear in the near future. 


3.4 - Other applications 


Às more information is becoming available on the 
chemistry, electrochemistry and engineering aspects 
concerning the behaviour of organic polymers, such 
as polyacetylene, more potential applications can be 
envisaged for these materials [47]. 

Apart from its conducting properties, important for 
the electrical and electronics industries, there is now 
a great deal of interest in improving the mechanical 
properties of polyacetylene. Therefore, one goal of 
recent research has been to find polyacetylene-alloys 
or composites having the electronic properties of 
doped polyacetylene and higher mechanical and 
optical properties [48, 49). These composites provide 
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protection for buildings, vehicles, etc., because they 
contain radar-absorbing particles, materials which 
provid IR camouflage at wavelengths = 14 u,and 
corrosion protection agents [50]. 

Recent studies have shown that (CH), acts as a 
“catalyst electrode”. It was shown that 05, H50s5 and 
HCl04 may be electrocatalytically reduced in 
aqueous electrolytes by (CH), [51], so new areas of 
applications (eg. fuel-cell type electrodes, 
electrochemical synthesis) exploiting this 
phenomenon, are appealing for attention. 
Electrically conductive porous synthetic polymeric 
compositions, comprising polyacetylene catalysts, 
have been obtained by U.S.A. workers [52], being 
useful for various membrane processes, namely 
electrodialysis. Other applications of polyacetylene 
include gas sensors [53], solid electrolytic capacitors 
[54], electroconductive inks [55], membranes for the 
separation of optically active compounds, selective 
chemical transducers for application in 
pharmacology and organic transistors [47]. 

In summary, the chemistry and physics of 
polyacetylene will open an extensive area for further 
theoretical and experimental research, not only of 
fundamental scientific interest but of potencial 
technological applications. 
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